A structural model for the assembly of the 30S subunit of the ribosome

Structural intermediates and folding events in the early assembly of the ribosomal small subunit

Using all-atom explicit solvent molecular dynamics (MD) simulations, we investigated the early structural intermediates of the 5' domain of the 16S rRNA in Escherichia coli upon the removal of the primary binding r-proteins S4, S17, and S20 and the secondary binding r-protein S16. Removal of each r-protein corresponded to the disappearance of subdomains with correlated dynamics. Correlation-based network analysis of the MD trajectories of the naked rRNA showed that the different subdomains are connected via multiple pathways with high betweenness. These pathways cross at the internal loop of helix 17 (h17) in the five-way junction (5WJ). The structure of the internal loop is disrupted by the binding of S17 and rescued by the addition of S16, suggesting an important function of the secondary binding protein in biasing the rRNA folding landscape toward the native basin. Using structure-based Gō simulations, we investigated the folding barriers of the lower four-way junction (4WJ) with h6, which is the primary binding site of S20 and the first to be transcribed. The folding of the 4WJ is consistent with the protection patterns observed in hydroxyl radical footprinting. Results from the all-atom simulations show that the fluctuations in the 5WJ are independent of the fluctuations in the 4WJ, suggesting that the subdomains fold independently and are stabilized by primary r-proteins.

RNA 3D structure prediction by using a coarse-grained model and experimental data

RNAs form complex secondary and three-dimensional structures, and their biological functions highly rely on their structures and dynamics. Here we developed a general coarse-grained framework for RNA 3D structure prediction. A new, hybrid coarse-grained model that explicitly describes the electrostatics and hydrogen-bond interactions has been constructed based on experimental structural statistics. With the simulated annealing simulation protocol, several RNAs of less than 30-nt were folded to within 4.0 Å of the native structures. In addition, with limited restraints on Watson-Crick basepairing based on the data from NMR spectroscopy and small-angle X-ray scattering (SAXS) information, the current model was able to characterize the complex tertiary structures of large size RNAs, such as 5S ribosome and U2/U6 snRNA. We also demonstrated that the pseudoknot structure was better captured when the coordinating Mg(2+) cations and limited basepairing restraints were included. The accuracy of our model has been compared favorably with other RNA structure prediction methods presented in the previous study of RNA-Puzzles. Therefore the coarse-grained model presented here offers a unique approach for accurate prediction and modeling of RNA structures.

A two-stage mechanism of viral RNA compaction revealed by single molecule fluorescence

Long RNAs often exist as multiple conformers in equilibrium. For the genomes of single-stranded RNA viruses, one of these conformers must include a compacted state allowing the RNA to be confined within the virion. We have used single molecule fluorescence correlation spectroscopy to monitor the conformations of viral genomes and sub-fragments in the absence and presence of coat proteins. Cognate RNA-coat protein interactions in two model viruses cause a rapid collapse in the hydrodynamic radii of their respective RNAs. This is caused by protein binding at multiple sites on the RNA that facilitate additional protein-protein contacts. The collapsed species recruit further coat proteins to complete capsid assembly with great efficiency and fidelity. The specificity in RNA-coat protein interactions seen at single-molecule concentrations reflects the packaging selectivity seen for such viruses in vivo. This contrasts with many in vitro reassembly measurements performed at much higher concentrations. RNA compaction by coat protein or polycation binding are distinct processes, implying that defined RNA-coat protein contacts are required for assembly.

Assembly of the five-way junction in the ribosomal small subunit using hybrid MD-Gō simulations

Assembly of the bacterial ribosomal small subunit (SSU) begins with the folding of the five-way junction upon interaction with the primary binding protein S4. This complex contains the largest contiguous molecular signature, which is a conserved feature in all bacterial 16S rRNAs. In a previous study, we used all-atom molecular dynamics simulations to demonstrate that the co-evolving signature in the N-terminus of S4 is intrinsically disordered and capable of accelerating the binding process through a fly casting mechanism. In this paper, comparisons between the all-atom MD simulations and FRET experiments identify multiple metastable conformations of the naked five-way junction without the presence of S4. Furthermore, we capture the simultaneous folding and binding of the five-way junction and r-protein S4 by use of a structure-based Gō potential implemented within the framework of the all-atom molecular dynamics CHARMM force field. Different folding pathways are observed for the refolding of the five-way junction upon partial binding of S4. Our simulations illustrate the complex nature of RNA folding in the presence of a protein binding partner and provide insight into the role of population shift and the induced fit mechanisms in the protein:RNA folding and binding process.

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.

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. Additionally, some antibacterial drugs are believed to target the bacterial ribosome during its construction. Hence, ribosomal assembly has been an active research topic for many years because understanding the assembly mechanisms can provide insight into protein/RNA recognitions important in many other cellular processes, as well as optimize the development of antibacterial therapeutics. Experimental studies thus far have provided still limited understanding about the assembly process. To further understand the assembly process, we have computationally studied the dynamic properties that r-proteins exhibit during assembly and the relationship between dynamics, physical properties, and binding propensity. We observe significant charged interactions between r-proteins and rRNA. We also detect a strong correlation between contact residues and their dynamic mobilities. Protein residues contacting with rRNA are observed to be more mobile in comparison with other residues. We also relate the location of the r-protein in the fully assembled ribosome to its susceptibility for large conformational changes prior to binding.

Coarse-grained models to study dynamics of nanoscale biomolecules and their applications to the ribosome

Biopolymers are of dynamic nature and undergo functional motions spanning a large spectrum of timescales. To study the internal dynamics of nano-sized molecular complexes that exceed hundred thousands of atoms with atomic detail is computationally inefficient. Therefore, to achieve both the spatial and temporal scales of biological interest coarse-grained models of macromolecules are often used. By uniting groups of atoms into single interacting centers one decreases the resolution of the system and gets rid of the irrelevant degrees of freedom. This simplification, even though it requires parameterization, makes the studies of biomolecular dynamics computationally tractable and allows us to reach beyond the microsecond time frame. Here, I review the coarse-grained models of macromolecules composed of proteins and nucleic acids. I give examples of one-bead models that were developed to investigate the internal dynamics and focus on their applications to the ribosome--the nanoscale protein synthesis machine.

Coarse-grained model for simulation of RNA three-dimensional structures

The accurate prediction of an RNA's three-dimensional structure from its "primary structure" will have a tremendous influence on the experimental design and its interpretation and ultimately our understanding of the many functions of RNA. This paper presents a general coarse-grained (CG) potential for modeling RNA 3-D structures. Each nucleotide is represented by five pseudo atoms, two for the backbone (one for the phosphate and another for the sugar) and three for the base to represent base-stacking interactions. The CG potential has been parametrized from statistical analysis of 688 RNA experimental structures. Molecular dynamic simulations of 15 RNA molecules with the length of 12-27 nucleotides have been performed using the CG potential, with performance comparable to that from all-atom simulations. For ~75% of systems tested, simulated annealing led to native-like structures at least once out of multiple repeated runs. Furthermore, with weak distance restraints based on the knowledge of three to five canonical Watson-Crick pairs, all 15 RNAs tested are successfully folded to within 6.5 Å of native structures using the CG potential and simulated annealing. The results reveal that with a limited secondary structure model the current CG potential can reliably predict the 3-D structures for small RNA molecules. We also explored an all-atom force field to construct atomic structures from the CG simulations.

Simulating activity of the bacterial ribosome

Computational modeling studies that investigate activity of the bacterial ribosome were reviewed. Computational approaches became possible with the availability of three-dimensional atomic resolution structures of the ribosomal subunits. However, due to the enormous size of the system, theoretical efforts to study the ribosome are few and challenging. For example, to extend the simulation timescales to biologically relevant ones, often, reduced models that require tedious parameterizations need to be applied. To that end, modeling of the ribosome focused on its internal dynamics, electrostatic properties, inhibition by antibiotics, polypeptide folding in the ribosome tunnel and assembly mechanisms driving the formation of the small ribosomal subunit.

On gene dosage balance in protein complexes: a comment on Semple JI, Vavouri T, Lehner B. A simple principle concerning the robustness of protein complex activity to changes in gene expression

A comment on Semple JI, Vavouri T, Lehner B. A simple principle concerning the robustness of protein complex activity to changes in gene expression. BMC Syst Biol. 2008;2:1

Effects of protein subunits removal on the computed motions of partial 30S structures of the ribosome

The Anisotropic Network Model (ANM) is used to study motions of the 30S small ribosomal subunit. The effect of the absence of certain subunits on the motions of the remaining partial structures was investigated by removing one protein, pairs of proteins and selected sets of proteins at a time. Our results show that the removal of some proteins doesn't change the large-scale dynamics of the partial structures, but the removal of certain subunits does cause significant changes in motion of the remaining structure, and these changes can be reverted by the removal of other subunits, which indicates interdependence between motions of various parts of the 30S ribosomal structure. We further found that the subunits showing such interdependence have strong positive correlation of their motions, which indicates that these subunits function as a unit block in the 30S small ribosomal subunit. Dynamically interdependent subunit pairs identified in this paper are consistent with previous experimental observations that suggested dimerization of those subunits.

Biophysical studies of bacterial ribosome assembly

The assembly of the bacterial ribosome involves the association of over 50 proteins to 3 large RNA molecules, and it represents a major metabolic activity for rapidly growing bacteria. The availability of atomic structures of the ribosome and the application of biochemical and biophysical methods have led to rapid progress in understanding the mechanistic details of ribosome assembly. The basic steps required to assemble a ribosome are outlined, and the contributions of mass spectrometry, computational methods, and RNA-folding studies in understanding these steps are detailed. This complex process takes place with both sequential and parallel processing that is coordinated to ensure efficient and complete assembly of ribosomes to meet the demands of cell growth.

Ribosome biogenesis and the translation process in Escherichia coli

Translation, the decoding of mRNA into protein, is the third and final element of the central dogma. The ribosome, a nucleoprotein particle, is responsible and essential for this process. The bacterial ribosome consists of three rRNA molecules and approximately 55 proteins, components that are put together in an intricate and tightly regulated way. When finally matured, the quality of the particle, as well as the amount of active ribosomes, must be checked. The focus of this review is ribosome biogenesis in Escherichia coli and its cross-talk with the ongoing protein synthesis. We discuss how the ribosomal components are produced and how their synthesis is regulated according to growth rate and the nutritional contents of the medium. We also present the many accessory factors important for the correct assembly process, the list of which has grown substantially during the last few years, even though the precise mechanisms and roles of most of the proteins are not understood.

Dependency map of proteins in the small ribosomal subunit

The assembly of the ribosome has recently become an interesting target for antibiotics in several bacteria. In this work, we extended an analytical procedure to determine native state fluctuations and contact breaking to investigate the protein stability dependence in the 30S small ribosomal subunit of Thermus thermophilus. We determined the causal influence of the presence and absence of proteins in the 30S complex on the binding free energies of other proteins. The predicted dependencies are in overall agreement with the experimentally determined assembly map for another organism, Escherichia coli. We found that the causal influences result from two distinct mechanisms: one is pure internal energy change, the other originates from the entropy change. We discuss the implications on how to target the ribosomal assembly most effectively by suggesting six proteins as targets for mutations or other hindering of their binding. Our results show that by blocking one out of this set of proteins, the association of other proteins is eventually reduced, thus reducing the translation efficiency even more. We could additionally determine the binding dependency of THX--a peptide not present in the ribosome of E. coli--and suggest its assembly path.

An assembly landscape for the 30S ribosomal subunit

Self-assembling macromolecular machines drive fundamental cellular processes, including transcription, messenger RNA processing, translation, DNA replication and cellular transport. The ribosome, which carries out protein synthesis, is one such machine, and the 30S subunit of the bacterial ribosome is the preeminent model system for biophysical analysis of large RNA-protein complexes. Our understanding of 30S assembly is incomplete, owing to the challenges of monitoring the association of many components simultaneously. Here we have developed a method involving pulse-chase monitored by quantitative mass spectrometry (PC/QMS) to follow the assembly of the 20 ribosomal proteins with 16S ribosomal RNA during formation of the functional particle. These data represent a detailed and quantitative kinetic characterization of the assembly of a large multicomponent macromolecular complex. By measuring the protein binding rates at a range of temperatures, we find that local transformations throughout the assembling subunit have similar but distinct activation energies. Thus, the prevailing view of 30S assembly as a pathway proceeding through a global rate-limiting conformational change must give way to one in which the assembly of the complex traverses a landscape dotted with various local conformational transitions.

Protein-independent Folding Pathway of the 16S rRNA 5′ Domain

Evolution of the ribosome from an RNA catalyst suggests that the intrinsic folding pathway of the rRNA dictates the hierarchy of ribosome assembly. To address this possibility, we probed the tertiary folding pathway of the 5' domain of the Escherichia coli 16S rRNA at 20 ms intervals using X-ray-dependent hydroxyl radical footprinting. Comparison with crystallographic structures and footprinting reactions on native 30S ribosomes showed that the RNA formed all of the predicted tertiary interactions in the absence of proteins. In 20 mM MgCl2, many tertiary interactions appeared within 20 ms. By contrast, interactions between H6, H15 and H17 near the spur of the 30S ribosome evolved over several minutes, likely due to mispairing of a central helix junction. The kinetic folding pathway of the RNA corresponded to the expected order of protein binding, suggesting that the RNA folding pathway forms the basis for early steps of ribosome assembly.

Exploring Assembly Energetics of the 30S Ribosomal Subunit Using an Implicit Solvent Approach

To explore the relationship between the assembly of the 30S ribosomal subunit and interactions among the constituent components, 16S RNA and proteins, relative binding free energies of the T. thermophilus 30S proteins to the 16S RNA were studied based on an implicit solvent model of electrostatic, nonpolar, and entropic contributions. The late binding proteins in our assembly map were found not to bind to the naked 16S RNA. The 5' domain early kinetic class proteins, on average, carry the highest positive charge, get buried the most upon binding to 16S RNA, and show the most favorable binding. Some proteins (S10/S14, S6/S18, S13/S19) have more stabilizing interactions while binding as dimers. Our computed assembly map resembles that of E. coli; however, the central domain path is more similar to that of A. aeolicus, a hyperthermophilic bacteria.

Nonbridging phosphate oxygens in 16S rRNA important for 30S subunit assembly and association with the 50S ribosomal subunit

Ribosomes are composed of RNA and protein molecules that associate together to form a supramolecular machine responsible for protein biosynthesis. Detailed information about the structure of the ribosome has come from the recent X-ray crystal structures of the ribosome and the ribosomal subunits. However, the molecular interactions between the rRNAs and the r-proteins that occur during the intermediate steps of ribosome assembly are poorly understood. Here we describe a modification-interference approach to identify nonbridging phosphate oxygens within 16S rRNA that are important for the in vitro assembly of the Escherichia coli 30S small ribosomal subunit and for its association with the 50S large ribosomal subunit. The 30S small subunit was reconstituted from phosphorothioate-substituted 16S rRNA and small subunit proteins. Active 30S subunits were selected by their ability to bind to the 50S large subunit and form 70S ribosomes. Analysis of the selected population shows that phosphate oxygens at specific positions in the 16S rRNA are important for either subunit assembly or for binding to the 50S subunit. The X-ray crystallographic structures of the 30S subunit suggest that some of these phosphate oxygens participate in r-protein binding, coordination of metal ions, or for the formation of intersubunit bridges in the mature 30S subunit. Interestingly, however, several of the phosphate oxygens identified in this study do not participate in any interaction in the mature 30S subunit, suggesting that they play a role in the early steps of the 30S subunit assembly.

Hinge-like motions in RNA kink-turns: the role of the second a-minor motif and nominally unpaired bases

Kink-turn (K-turn) motifs are asymmetric internal loops found at conserved positions in diverse RNAs, with sharp bends in phosphodiester backbones producing V-shaped structures. Explicit-solvent molecular dynamics simulations were carried out for three K-turns from 23S rRNA, i.e., Kt-38 located at the base of the A-site finger, Kt-42 located at the base of the L7/L12 stalk, and Kt-58 located in domain III, and for the K-turn of human U4 snRNA. The simulations reveal hinge-like K-turn motions on the nanosecond timescale. The first conserved A-minor interaction between the K-turn stems is entirely stable in all simulations. The angle between the helical arms of Kt-38 and Kt-42 is regulated by local variations of the second A-minor (type I) interaction between the stems. Its variability ranges from closed geometries to open ones stabilized by insertion of long-residency waters between adenine and cytosine. The simulated A-minor geometries fully agree with x-ray data. Kt-58 and Kt-U4 exhibit similar elbow-like motions caused by conformational change of the adenosine from the nominally unpaired region. Despite the observed substantial dynamics of K-turns, key tertiary interactions are stable and no sign of unfolding is seen. We suggest that some K-turns are flexible elements mediating large-scale ribosomal motions during the protein synthesis cycle.

A determination of the identity elements in yeast 18 S ribosomal RNA for the recognition of ribosomal protein YS11: the role of the kink-turn motif in helix 11

A description of the site of interaction of YS11, the yeast homolog of eubacterial S17, with 18 S rRNA was obtained by assessing the binding of the ribosomal protein, in a filter retention assay, to oligoribonucleotides that reproduce regions of 18 S rRNA. YS11 binds predominantly to domain I; the Kd value is 113 nM. The dimensions of the YS11 binding site were refined, guided by chemical protection data and by the atomic structure of the Thermus thermophilus 30 S subunit, which has the S17 recognition site in 16 S rRNA. An oligoribonucleotide that mimics helix 11, a phylogenetically conserved region in domain I, binds YS11 with a Kd value of 230 nM; a second oligoribonucleotide that contains only the kink-turn motif in helix 11 binds YS11 with a Kd value of 528 nM. Thus, helix 11 has most of the nucleotides required for the recognition of YS11. To identify those nucleotides a set of 27 transversion mutations in H11 was constructed and their contribution to the binding of YS11 determined. Mutations of nine nucleotides (U313, C314, A316, G337, C338, G347, U348, U350, and C351) increased the Kd value for YS11 binding by at least eightfold; G325U and U349A mutations increased the Kd value fivefold. Eight of the 11 mutations are in the kink-turn in H11, confirming the critical importance of the motif for YS11 recognition. The other three nucleotides are in the lower stem and the terminal loop of H11, which makes a lesser, but still important, contribution to YS11 binding. The identity elements for YS11 recognition are: A316, G325, G337, G347, U348, U349, U350, and C351. The effect of the other nucleotides that decrease binding is probably indirect, presumably they affect the conformation of the binding site but do not have contacts to YS11 amino acid residues. The eight identity element nucleotides are in regions of H11 that deviate from A-form geometry and the contacts are predominantly, if not exclusively, to backbone phosphate and sugar oxygen atoms, indicating that YS11 recognizes the shape of the rRNA binding site rather than reading the sequence of nucleotides.

The YrdC protein--a putative ribosome maturation factor

Release factor one (RF1) terminates protein synthesis in response to stop codons UAG and UAA. A mutant allele of RF1 causes temperature sensitive growth at 42 degrees C. We have earlier described the isolation of a suppressor of the temperature sensitive phenotype. The suppressor mutation is a small deletion in the open reading frame yrdC, and we have shown that the DeltayrdC mutation leads to immature 30S subunits and, as a consequence, to fewer translating ribosomes. YrdC is a small conserved protein with a dsRNA-binding surface. Here, we have characterized the YrdC protein. We show that the deletion leads to no production of functional protein, and we have indications that the YrdC protein might be essential in a wild type background. The protein is needed for the maturation of 16S rRNA, even though it does not interact tightly with either of the ribosomal subunits, or the 70S particles. The less effective maturation of rRNA affects the ribosomal feedback control, leading to an increase in expression from P1rrnB. We suggest that the function of the YrdC protein is to keep an rRNA structure needed for proper processing of 16S rRNA, especially at lower temperatures. This activity may require other factor(s). We suggest the gene be renamed rimN, and the mutant allele rimN141.

Automated identification of RNA conformational motifs: theory and application to the HM LSU 23S rRNA

We develop novel methods for recognizing and cataloging conformational states of RNA, and for discovering statistical rules governing those states. We focus on the conformation of the large ribosomal subunit from Haloarcula marismortui. The two approaches described here involve torsion matching and binning. Torsion matching is a pattern-recognition code which finds structural repetitions. Binning is a classification technique based on distributional models of the data. In comparing the results of the two methods we have tested the hypothesis that the conformation of a very large complex RNA molecule can be described accurately by a limited number of discrete conformational states. We identify and eliminate extraneous and redundant information without losing accuracy. We conclude, as expected, that four of the torsion angles contain the overwhelming bulk of the structural information. That information is not significantly compromised by binning the continuous torsional information into a limited number of discrete values. The correspondence between torsion matching and binning is 99% (per residue). Binning, however, does have several advantages. In particular, we demonstrate that the conformation of a large complex RNA molecule can be represented by a small alphabet. In addition, the binning method lends itself to a natural graphical representation using trees.

Authentic precursors to ribosomal subunits accumulate in Escherichia coli in the absence of functional DnaK chaperone

Escherichia coli dnaK-ts mutants are defective in the late stages of ribosome biogenesis at high temperature. Here, we show that the 21S, 32S and 45S ribosomal particles that accumulate in the dnaK756-ts mutant at 44 degrees C contain unprocessed forms of their 16S and 23S rRNAs (partially processed in the case of 45S particles). Their 5S rRNA stoichiometry and ribosomal protein composition are typical of the genuine ribosomal precursors found in a wild-type (dnaK+) strain. Despite the lack of a functional DnaK, a very slow maturation of these 21S, 32S and 45S particles to structurally and functionally normal 30S and 50S ribosomal subunits still occurs at high temperature. This conversion is accompanied by the processing of p16S and p23S rRNAs to their mature forms. We conclude that: (i) 21S, 32S and 45S particles are not dead-end particles, but true precursors to active ribosomes (21S particles are converted to 30S subunits, and 32S and 45S to 50S subunits); (ii) DnaK is not absolutely necessary for ribosome biogenesis, but accelerates the late steps of this process considerably at high temperature; and (iii) 23S rRNA processing depends on the stage reached in the stepwise assembly of the 50S subunit, not directly on DnaK.

Ribosomal dynamics inferred from variations in experimental measurements

The crystal structures of the ribosome reveal remarkable complexity and provide a starting set of snapshots with which to understand the dynamics of translation. To augment the static crystallographic models with dynamic information present in crosslink, footprint, and cleavage data, we examined 2691 proximity measurements and focused on the subset that was apparently incompatible with >40 published crystal structures. The measurements from this subset generally involve regions of the structure that are functionally conserved and structurally flexible. Local movements in the crystallographic states of the ribosome that would satisfy biochemical proximity measurements show coherent patterns suggesting alternative conformations of the ribosome. Three different types of data obtained for the two subunits display similar "mismatching" patterns, suggesting that the signals are robust and real. In particular, there is an indication of coherent motion in the decoding region within the 30S subunit and central protuberance and surrounding areas of the 50S subunit. Directions of rearrangements fluctuate around the proposed path of tRNA translocation and the plane parallel to the interface of the two subunits. Our results demonstrate that systematic combination and analysis of noisy, apparently incompatible data sources can provide biologically useful signals about structural dynamics.