The Native 3D Organization of Bacterial Polysomes

SummaryRecent advances have led to insights into the structure of the bacterial ribosome, but little is known about the 3D organization of ribosomes in the context of translating polysomes. We employed cryoelectron tomography and a template-matching approach to map 70S ribosomes in vitrified bacterial translation extracts and in lysates of active E. coli spheroplasts. In these preparations, polysomal arrangements were observed in which neighboring ribosomes are densely packed and exhibit preferred orientations. Analysis of characteristic examples of polysomes reveals a staggered or pseudohelical organization of ribosomes along the mRNA trace, with the transcript being sequestered on the inside, the tRNA entrance sites being accessible, and the polypeptide exit sites facing the cytosol. Modeling of elongating nascent polypeptide chains suggests that this arrangement maximizes the distance between nascent chains on adjacent ribosomes, thereby reducing the probability of intermolecular interactions that would give rise to aggregation and limit productive folding

Diffusion, Crowding & Protein Stability in a Dynamic Molecular Model of the Bacterial Cytoplasm

A longstanding question in molecular biology is the extent to which the behavior of macromolecules observed in vitro accurately reflects their behavior in vivo. A number of sophisticated experimental techniques now allow the behavior of individual types of macromolecule to be studied directly in vivo; none, however, allow a wide range of molecule types to be observed simultaneously. In order to tackle this issue we have adopted a computational perspective, and, having selected the model prokaryote Escherichia coli as a test system, have assembled an atomically detailed model of its cytoplasmic environment that includes 50 of the most abundant types of macromolecules at experimentally measured concentrations. Brownian dynamics (BD) simulations of the cytoplasm model have been calibrated to reproduce the translational diffusion coefficients of Green Fluorescent Protein (GFP) observed in vivo, and “snapshots” of the simulation trajectories have been used to compute the cytoplasm's effects on the thermodynamics of protein folding, association and aggregation events. The simulation model successfully describes the relative thermodynamic stabilities of proteins measured in E. coli, and shows that effects additional to the commonly cited “crowding” effect must be included in attempts to understand macromolecular behavior in vivo

Cotranslational protein assembly imposes evolutionary constraints on homomeric proteins

Cotranslational protein folding can facilitate rapid formation of functional structures. However, it might also cause premature assembly of protein complexes, if two interacting nascent chains are in close proximity. By analyzing known protein structures, we show that homomeric protein contacts are enriched towards the C-termini of polypeptide chains across diverse proteomes. We hypothesize that this is the result of evolutionary constraints for folding to occur prior to assembly. Using high-throughput imaging of protein homomers in vivo in E. coli and engineered protein constructs with N- and C-terminal oligomerization domains, we show that, indeed, proteins with C-terminal homomeric interface residues consistently assemble more efficiently than those with N-terminal interface residues. Using in vivo, in vitro and in silico experiments, we identify features that govern successful assembly of homomers, which have implications for protein design and expression optimization

Molecular simulations of cotranslational protein folding: fragment stabilities, folding cooperativity, and trapping in the ribosome.

Although molecular simulation methods have yielded valuable insights into mechanistic aspects of protein refolding in vitro, they have up to now not been used to model the folding of proteins as they are actually synthesized by the ribosome. To address this issue, we report here simulation studies of three model proteins: chymotrypsin inhibitor 2 (CI2), barnase, and Semliki forest virus protein (SFVP), and directly compare their folding during ribosome-mediated synthesis with their refolding from random, denatured conformations. To calibrate the methodology, simulations are first compared with in vitro data on the folding stabilities of N-terminal fragments of CI2 and barnase; the simulations reproduce the fact that both the stability and thermal folding cooperativity increase as fragments increase in length. Coupled simulations of synthesis and folding for the same two proteins are then described, showing that both fold essentially post-translationally, with mechanisms effectively identical to those for refolding. In both cases, confinement of the nascent polypeptide chain within the ribosome tunnel does not appear to promote significant formation of native structure during synthesis; there are however clear indications that the formation of structure within the nascent chain is sensitive to location within the ribosome tunnel, being subject to both gain and loss as the chain lengthens. Interestingly, simulations in which CI2 is artificially stabilized show a pronounced tendency to become trapped within the tunnel in partially folded conformations: non-cooperative folding, therefore, appears in the simulations to exert a detrimental effect on the rate at which fully folded conformations are formed. Finally, simulations of the two-domain protease module of SFVP, which experimentally folds cotranslationally, indicate that for multi-domain proteins, ribosome-mediated folding may follow different pathways from those taken during refolding. Taken together, these studies provide a first step toward developing more realistic methods for simulating protein folding as it occurs in vivo

Atomistic simulations of competition between substrates binding to an enzyme.

Although the idea that electrostatic potentials generated by enzymes can guide substrates to active sites is well established, it is not always appreciated that the same potentials can also promote the binding of molecules other than the intended substrate, with the result that such enzymes might be sensitive to the presence of competing molecules. To provide a novel means of studying such "electrostatic competition" effects, computer simulation methodology has been developed to allow the diffusion and association of many solute molecules around a single enzyme to be simulated. To demonstrate the power of the methodology, simulations have been conducted on an artificial fusion protein of citrate synthase (CS) and malate dehydrogenase (MDH) to assess the chances of oxaloacetate being channeled between the MDH and CS active sites. The simulations demonstrate that the probability of channeling is strongly dependent on the concentration of the initial substrate (malate) in the solution. In fact, the high concentrations of malate used in experiments appear high enough to abolish any channeling of oxaloacetate. The simulations provide a resolution of a serious discrepancy between previous simulations and experiments and raise important questions relating to the observability of electrostatically mediated substrate channeling in vitro and in vivo

Parameter Sensitivity of Coupled Synthesis–Folding Simulations of CI2

<div><p>(A) Mean values of Q during coupled synthesis–folding of CI2 comparing the effects of using different timescales for amino acid synthesis (τ_syn).</p><p>(B) Same as (A), but comparing the effects of assigning larger radii to ribosome pseudo-atoms.</p><p>(C) Same as (A), but comparing the effects of including energetically favorable interactions between the nascent protein and the ribosome.</p><p>(D) Distribution of mean Q_res values obtained during the synthesis of the final amino acid, comparing the effects of assigning larger radii to ribosome pseudo-atoms; error bars indicate the standard deviation of values obtained from 30 independent trajectories.</p><p>(E) Same as (D), but comparing the effects of including energetically favorable interactions between the nascent protein and the ribosome.</p></div

Consequences of Excessively Stabilizing CI2

<div><p>(A) Histograms of folding times obtained from refolding and de novo folding trajectories of CI2 with ɛ = 0.80 kcal/mol.</p><p>(B) Snapshot of a “trapped” trajectory.</p><p>(C) Plots of Q versus time for ten previously trapped trajectories when restarted without the ribosome.</p><p>(D) Same as (C), but restarted with the ribosome and ɛ = 0.60 kcal/mol.</p></div

A Typical Coupled Synthesis–Folding Simulation of SFVP

<p>Snapshots from a typical de novo folding simulation of SFVP in the presence of the ribosome (grey surface). The nascent chain is colored from blue (N-terminus) to red (C-terminus), with its length (in amino acids) indicated in each panel. This figure was prepared with PyMOL[<a href="http://www.ploscompbiol.org/article/info:doi/10.1371/journal.pcbi.0020098#pcbi-0020098-b095" target="_blank">95</a>].</p