Ribosome motions modulate electrostatic properties

Defining the EM-signature of successful cell-transfection

In this report, we describe the architecture of Lipofectamine 2000 and 3000 transfection- reagents, as they appear inside of transfected cells, using classical transmission electron microscopy (EM). We also demonstrate that they provoke consistent structural changes after they have entered cells, changes that not only provide new insights into the mechanism of action of these particular transfection-reagents, but also provide a convenient and robust method for identifying by EM which cells in any culture have been successfully transfected. This also provides clues to the mechanism(s) of their toxic effects, when they are applied in excess. We demonstrate that after being bulk-endocytosed by cells, the cationic spheroids of Lipofectamine remain intact throughout the entire time of culturing, but escape from their endosomes and penetrate directly into the cytoplasm of the cell. In so doing, they provoke a stereotypical recruitment and rearrangement of endoplasmic reticulum (ER), and they ultimately end up escaping into the cytoplasm and forming unique 'inclusion-bodies.' Once free in the cytoplasm, they also invariably develop dense and uniform coatings of cytoplasmic ribosomes on their surfaces, and finally, they become surrounded by 'annulate' lamellae' of the ER. In the end, these annulate-lamellar enclosures become the ultrastructural 'signatures' of these inclusion-bodies, and serve to positively and definitively identify all cells that have been effectively transfected. Importantly, these new EM-observations define several new and unique properties of these classical Lipofectamines, and allow them to be discriminated from other lipoidal or particulate transfection-reagents, which we find do not physically break out of endosomes or end up in inclusion bodies, and in fact, provoke absolutely none of these 'signature' cytoplasmic reactions.

Electrostatic Ratchet for Successive Peptide Synthesis in Nonribosomal Molecular Machine RimK

A nonribosomal peptide-synthesizing molecular machine, RimK, adds l-glutamic acids to the C-terminus of ribosomal protein S6 (RpsF) in vivo and synthesizes poly-α-glutamates in vitro. However, the mechanism of the successive glutamate addition, which is fueled by ATP, remains unclear. Here, we investigate the successive peptide-synthesizing mechanism of RimK via the molecular dynamics (MD) simulation of glutamate binding. We first show that RimK adopts three stable structural states with respect to the ATP-binding loop and the triphosphate chain of the bound ATP. We then show that a glutamate in solution preferentially binds to a positively charged belt-like region of RimK and the bound glutamate exhibits Brownian motion along the belt. The binding-energy landscape shows that the open-to-closed transition of the ATP-binding loop and the bent-to-straight transition of the triphosphate chain of ATP can function as an electrostatic ratchet that guides the bound glutamate to the active site. We then show the binding site of the second glutamate, which allows us to infer the ligation mechanism. Consistent with MD results, the crystal structure of RimK we obtained in the presence of RpsF presents an electron density that is presumed to correspond to the C-terminus of RpsF. We finally propose a mechanism for the successive peptide synthesis by RimK and discuss its similarity to other molecular machines.

The significance of bioelectricity on all levels of organization of an organism. Part 1: From the subcellular level to cells

Bioelectricity plays an essential role in the structural and functional organization of biological organisms. In this first article of our three-part series, we summarize the importance of bioelectricity for the basic structural level of biological organization, i.e. from the subcellular level (charges, ion channels, molecules and cell organelles) to cells.

The energetics of subunit rotation in the ribosome

Protein synthesis in the cell is controlled by an elaborate sequence of conformational rearrangements in the ribosome. The composition of a ribosome varies by species, though they typically contain ∼ 50-100 RNA and protein molecules. While advances in structural techniques have revolutionized our understanding of long-lived conformational states, a vast range of transiently visited configurations can not be directly observed. In these cases, computational/simulation methods can be used to understand the mechanical properties of the ribosome. Insights from these approaches can then help guide next-generation experimental measurements. In this short review, we discuss theoretical strategies that have been deployed to quantitatively describe the energetics of collective rearrangements in the ribosome. We focus on efforts to probe large-scale subunit rotation events, which involve the coordinated displacement of large numbers of atoms (tens of thousands). These investigations are revealing how the molecular structure of the ribosome encodes the mechanical properties that control large-scale dynamics.

In silico analysis of bacterial translation factors reveal distinct translation event specific pI values

Background

Protein synthesis is a cellular process that takes place through the successive translation events within the ribosome by the event-specific protein factors, namely, initiation, elongation, release, and recycling factors. In this regard, we asked the question about how similar are those translation factors to each other from a wide variety of bacteria? Hence, we did a thorough in silico study of the translation factors from 495 bacterial sp., and 4262 amino acid sequences by theoretically measuring their pI and MW values that are two determining factors for distinguishing individual proteins in 2D gel electrophoresis in experimental procedures. Then we analyzed the output from various angles.

Results

Our study revealed the fact that it’s not all same, or all random, but there are distinct orders and the pI values of translation factors are translation event specific. We found that the translation initiation factors are mainly basic, whereas, elongation and release factors that interact with the inter-subunit space of the intact 70S ribosome during translation are strictly acidic across bacterial sp. These acidic elongation factors and release factors contain higher frequencies of glutamic acids. However, among all the translation factors, the translation initiation factor 2 (IF2) and ribosome recycling factor (RRF) showed variable pI values that are linked to the order of phylogeny.

Conclusions

From the results of our study, we conclude that among all the bacterial translation factors, elongation and release factors are more conserved in terms of their pI values in comparison to initiation and recycling factors. Acidic properties of these factors are independent of habitat, nature, and phylogeny of the bacterial species. Furthermore, irrespective of the different shapes, sizes, and functions of the elongation and release factors, possession of the strictly acidic pI values of these translation factors all over the domain Bacteria indicates that the acidic nature of these factors is a necessary criterion, perhaps to interact into the partially enclosed rRNA rich inter-subunit space of the translating 70S ribosome.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07472-x.

Charge-driven condensation of RNA and proteins suggests broad role of phase separation in cytoplasmic environments

Phase separation processes are increasingly being recognized as important organizing mechanisms of biological macromolecules in cellular environments. Well-established drivers of phase separation are multi-valency and intrinsic disorder. Here, we show that globular macromolecules may condense simply based on electrostatic complementarity. More specifically, phase separation of mixtures between RNA and positively charged proteins is described from a combination of multiscale computer simulations with microscopy and spectroscopy experiments. Phase diagrams were mapped out as a function of molecular concentrations in experiment and as a function of molecular size and temperature via simulations. The resulting condensates were found to retain at least some degree of internal dynamics varying as a function of the molecular composition. The results suggest a more general principle for phase separation that is based primarily on electrostatic complementarity without invoking polymer properties as in most previous studies. Simulation results furthermore suggest that such phase separation may occur widely in heterogenous cellular environment between nucleic acid and protein components.

Complex long-distance effects of mutations that confer linezolid resistance in the large ribosomal subunit

The emergence of multidrug-resistant pathogens will make current antibiotics ineffective. For linezolid, a member of the novel oxazolidinone class of antibiotics, 10 nucleotide mutations in the ribosome have been described conferring resistance. Hypotheses for how these mutations affect antibiotics binding have been derived based on comparative crystallographic studies. However, a detailed description at the atomistic level of how remote mutations exert long-distance effects has remained elusive. Here, we show that the G2032A-C2499A double mutation, located > 10 Å away from the antibiotic, confers linezolid resistance by a complex set of effects that percolate to the binding site. By molecular dynamics simulations and free energy calculations, we identify U2504 and C2452 as spearheads among binding site nucleotides that exert the most immediate effect on linezolid binding. Structural reorganizations within the ribosomal subunit due to the mutations are likely associated with mutually compensating changes in the effective energy. Furthermore, we suggest two main routes of information transfer from the mutation sites to U2504 and C2452. Between these, we observe cross-talk, which suggests that synergistic effects observed for the two mutations arise in an indirect manner. These results should be relevant for the development of oxazolidinone derivatives that are active against linezolid-resistant strains.

Using Molecular Simulation to Model High-Resolution Cryo-EM Reconstructions

An explosion of new data from high-resolution cryo-electron microscopy (cryo-EM) studies has produced a large number of data sets for many species of ribosomes in various functional states over the past few years. While many methods exist to produce structural models for lower resolution cryo-EM reconstructions, high-resolution reconstructions are often modeled using crystallographic techniques and extensive manual intervention. Here, we present an automated fitting technique for high-resolution cryo-EM data sets that produces all-atom models highly consistent with the EM density. Using a molecular dynamics approach, atomic positions are optimized with a potential that includes the cross-correlation coefficient between the structural model and the cryo-EM electron density, as well as a biasing potential preserving the stereochemistry and secondary structure of the biomolecule. Specifically, we use a hybrid structure-based/ab initio molecular dynamics potential to extend molecular dynamics fitting. In addition, we find that simulated annealing integration, as opposed to straightforward molecular dynamics integration, significantly improves performance. We obtain atomistic models of the human ribosome consistent with high-resolution cryo-EM reconstructions of the human ribosome. Automated methods such as these have the potential to produce atomistic models for a large number of ribosome complexes simultaneously that can be subsequently refined manually.

The effect of tRNA levels on decoding times of mRNA codons

The possible effect of transfer ribonucleic acid (tRNA) concentrations on codons decoding time is a fundamental biomedical research question; however, due to a large number of variables affecting this process and the non-direct relation between them, a conclusive answer to this question has eluded so far researchers in the field. In this study, we perform a novel analysis of the ribosome profiling data of four organisms which enables ranking the decoding times of different codons while filtering translational phenomena such as experimental biases, extreme ribosomal pauses and ribosome traffic jams. Based on this filtering, we show for the first time that there is a significant correlation between tRNA concentrations and the codons estimated decoding time both in prokaryotes and in eukaryotes in natural conditions (−0.38 to −0.66, all P values <0.006); in addition, we show that when considering tRNA concentrations, codons decoding times are not correlated with aminoacyl-tRNA levels. The reported results support the conjecture that translation efficiency is directly influenced by the tRNA levels in the cell. Thus, they should help to understand the evolution of synonymous aspects of coding sequences via the adaptation of their codons to the tRNA pool.

A comparative study of ribosomal proteins: linkage between amino acid distribution and ribosomal assembly

BACKGROUND

Assembly of the ribosome from its protein and RNA constituents must occur quickly and efficiently in order to synthesize the proteins necessary for all cellular activity. Since the early 1960's, certain characteristics of possible assembly pathways have been elucidated, yet the mechanisms that govern the precise recognition events remain unclear.We utilize a comparative analysis to investigate the amino acid composition of ribosomal proteins (r-proteins) with respect to their role in the assembly process. We compared small subunit (30S) r-protein sequences to those of other housekeeping proteins from 560 bacterial species and searched for correlations between r-protein amino acid content and factors such as assembly binding order, environmental growth temperature, protein size, and contact with ribosomal RNA (rRNA) in the 30S complex.

RESULTS

We find r-proteins have a significantly high percent of positive residues, which are highly represented at rRNA contact sites. An inverse correlation between the percent of positive residues and r-protein size was identified and is mainly due to the content of Lysine residues, rather than Arginine. Nearly all r-proteins carry a net positive charge, but no statistical correlation between the net charge and the binding order was detected. Thermophilic (high-temperature) r-proteins contain increased Arginine, Isoleucine, and Tyrosine, and decreased Serine and Threonine compared to mesophilic (lower-temperature), reflecting a known distinction between thermophiles and mesophiles, possibly to account for protein thermostability. However, this difference in amino acid content does not extend to rRNA contact sites, as the proportions of thermophilic and mesophilic contact residues are not significantly different.

CONCLUSIONS

Given the significantly higher level of positively charged residues in r-proteins and at contact sites, we conclude that ribosome assembly relies heavily on an electrostatic component of interaction. However, the binding order of r-proteins in assembly does not appear to depend on these electrostatics interactions. Additionally, because thermophiles and mesophiles exhibit significantly different amino acid compositions in their sequences but not in the identities of contact sites, we conclude that this electrostatic component of interaction is insensitive to temperature and is not the determining factor differentiating the temperature sensitivity of ribosome assembly.

Electrostatic effect of the ribosomal surface on nascent polypeptide dynamics

The crucial molecular events accompanying protein folding in the cell are still largely unexplored. As nascent polypeptides emerge from the ribosomal exit tunnel, they come in close proximity with the highly negatively charged ribosomal surface. How is the nascent polypeptide influenced by the ribosomal surface? We address this question via the intrinsically disordered protein PIR and a number of its variably charged mutants. Two different populations are identified: one is highly spatially biased, and the other is highly dynamic. The more negatively charged nascent polypeptides emerging from the ribosome are richer in the extremely dynamic population. Hence, nascent proteins with a net negative charge are less likely to interact with the ribosome. Surprisingly, the amplitude of the local motions of the highly dynamic population is much wider than that of disordered polypeptides under physiological conditions, implying that proximity to the ribosomal surface enhances the molecular flexibility of a subpopulation of the nascent protein, much like a denaturing agent would. This effect could be important for a proper structural channeling of the nascent protein and the prevention of cotranslational kinetic trapping. Interestingly, a significant population of the highly spatially biased nascent chain, probably interacting extensively with the ribosome, is present even for very negatively charged nascent proteins. This "sticking" effect likely serves to protect nascent proteins (e.g., from cotranslational aggregation). In all, our results highlight the influence of the ribosome in nascent protein dynamics and show that the ribosome's function in protein biogenesis extends well beyond catalysis of peptide bond formation.

Determinants of translation elongation speed and ribosomal profiling biases in mouse embryonic stem cells

Ribosomal profiling is a promising approach with increasing popularity for studying translation. This approach enables monitoring the ribosomal density along genes at a resolution of single nucleotides.

In this study, we focused on ribosomal density profiles of mouse embryonic stem cells. Our analysis suggests, for the first time, that even in mammals such as M. musculus the elongation speed is significantly and directly affected by determinants of the coding sequence such as: 1) the adaptation of codons to the tRNA pool; 2) the local mRNA folding of the coding sequence; 3) the local charge of amino acids encoded in the codon sequence. In addition, our analyses suggest that in general, the translation velocity of ribosomes is slower at the beginning of the coding sequence and tends to increase downstream.

Finally, a comparison of these data to the expected biophysical behavior of translation suggests that it suffers from some unknown biases. Specifically, the ribosomal flux measured on the experimental data increases along the coding sequence; however, according to any biophysical model of ribosomal movement lacking internal initiation sites, the flux is expected to remain constant or decrease. Thus, developing experimental and/or statistical methods for understanding, detecting and dealing with such biases is of high importance.

Gene translation is the process by which ribosomes translate mRNA molecules to proteins, a central process in all living organisms. Thus, understanding the biophysics of gene translation and the way its efficiency is encoded in the different features of the coding sequence has ramifications to every biomedical discipline. Recently, a new large-scale experimental approach named ‘ribosomal profiling’, has been developed for monitoring the ribosomal density at a resolution of single nucleotides. In this study, we analyzed ribosomal profiling data of mouse embryonic stem cells. These data enabled us to directly show that translation velocity is affected by the adaptation of codons to the tRNA pool, local mRNA folding of coding sequence, and local charge of the amino acids encoded in the coding sequence. In addition, our analyses suggest that ribosomal speed tends to be slower at the beginning of the coding sequence. Finally, we report possible biases in the ‘ribosomal profiling’ procedure that should be considered in future studies utilizing this method.

iAPBS: a programming interface to Adaptive Poisson-Boltzmann Solver (APBS)

The Adaptive Poisson-Boltzmann Solver (APBS) is a state-of-the-art suite for performing Poisson-Boltzmann electrostatic calculations on biomolecules. The iAPBS package provides a modular programmatic interface to the APBS library of electrostatic calculation routines. The iAPBS interface library can be linked with a FORTRAN or C/C++ program thus making all of the APBS functionality available from within the application. Several application modules for popular molecular dynamics simulation packages - Amber, NAMD and CHARMM are distributed with iAPBS allowing users of these packages to perform implicit solvent electrostatic calculations with APBS.

Unwrapping of nucleosomal DNA ends: a multiscale molecular dynamics study

To permit access to DNA-binding proteins involved in the control and expression of the genome, the nucleosome undergoes structural remodeling including unwrapping of nucleosomal DNA segments from the nucleosome core. Here we examine the mechanism of DNA dissociation from the nucleosome using microsecond timescale coarse-grained molecular dynamics simulations. The simulations exhibit short-lived, reversible DNA detachments from the nucleosome and long-lived DNA detachments not reversible on the timescale of the simulation. During the short-lived DNA detachments, 9 bp dissociate at one extremity of the nucleosome core and the H3 tail occupies the space freed by the detached DNA. The long-lived DNA detachments are characterized by structural rearrangements of the H3 tail including the formation of a turn-like structure at the base of the tail that sterically impedes the rewrapping of DNA on the nucleosome surface. Removal of the H3 tails causes the long-lived detachments to disappear. The physical consistency of the CG long-lived open state was verified by mapping a CG structure representative of this state back to atomic resolution and performing molecular dynamics as well as by comparing conformation-dependent free energies. Our results suggest that the H3 tail may stabilize the nucleosome in the open state during the initial stages of the nucleosome remodeling process.

Computational studies of molecular machines: the ribosome

The past decade has produced an avalanche of experimental data on the structure and dynamics of the ribosome. Groundbreaking studies in structural biology and kinetics have placed important constraints on ribosome structural dynamics. However, a gulf remains between static structures and time dependent data. In particular, X-ray crystallography and cryo-EM studies produce static models of the ribosome in various states, but lack dynamic information. Single molecule studies produce information on the rates of transitions between these states but do not have high-resolution spatial information. Computational studies have aided in bridging this gap by providing atomic resolution simulations of structural fluctuations and transitions between configurations.

On the origin of order in the genome organization of ssRNA viruses

Single-stranded RNA (ssRNA) viruses form a major class that includes important human, animal, and plant pathogens. While the principles underlying the structures of their protein capsids are generally well understood, much less is known about the organization of their encapsulated genomic RNAs. Cryo-electron microscopy and x-ray crystallography have revealed striking evidence of order in the packaged genomes of a number of ssRNA viruses. The physical determinants of such order, however, are largely unknown. We study here the relative effect of different energetic contributions, as well as the role of confinement, on the genome packaging of a representative ssRNA virus, the bacteriophage MS2, via a series of biomolecular simulations in which different energy terms are systematically switched off. We show that the bimodal radial density profile of the packaged genome is a consequence of RNA self-repulsion in confinement, suggesting that it should be similar for all ssRNA viruses with a comparable ratio of capsid size/genome length. In contrast, the detailed structure of the outer shell of the RNA density depends crucially on steric contributions from the capsid inner surface topography, implying that the various different polyhedral RNA cages observed in experiment are largely due to differences in the inner surface topography of the capsid.

Composite effects of gene determinants on the translation speed and density of ribosomes

Background

Translation is a central process of life, and its regulation is crucial for cell growth. In this article, focusing on two model organisms, Escherichia coli and Saccharomyces cerevisiae, we study how three major local features of a gene's coding sequence (its adaptation to the tRNA pool, its amino acid charge, and its mRNA folding energy) affect its translation elongation.

Results

We find that each of these three different features has a non-negligible distinct correlation with the speed of translation elongation. In addition, each of these features might contribute independently to slowing down ribosomal speed at the beginning of genes, which was suggested in previous studies to improve ribosomal allocation and the cost of translation, and to decrease ribosomal jamming. Remarkably, a model of ribosomal translation based on these three basic features highly correlated with the genomic profile of ribosomal density. The robustness to transcription errors in terms of the values of these features is higher at the beginnings of genes, suggesting that this region is important for translation.

Conclusions

The reported results support the conjecture that translation elongation speed is affected by the three coding sequence determinants mentioned above, and not only by adaptation to the tRNA pool; thus, evolution shapes all these determinants along the coding sequences and across genes to improve the organism's translation efficiency.

Brownian dynamics study of the association between the 70S ribosome and elongation factor G

Protein synthesis on the ribosome involves a number of external protein factors that bind at its functional sites. One key factor is the elongation factor G (EF-G) that facilitates the translocation of transfer RNAs between their binding sites, as well as advancement of the messenger RNA by one codon. The details of the EF-G/ribosome diffusional encounter and EF-G association pathway still remain unanswered. Here, we applied Brownian dynamics methodology to study bimolecular association in the bacterial EF-G/70S ribosome system. We estimated the EF-G association rate constants at 150 and 300 mM monovalent ionic strengths and obtained reasonable agreement with kinetic experiments. We have also elucidated the details of EF-G/ribosome association paths and found that positioning of the L11 protein of the large ribosomal subunit is likely crucial for EF-G entry to its binding site.

A Fast and Robust Poisson-Boltzmann Solver Based on Adaptive Cartesian Grids

An adaptive Cartesian grid (ACG) concept is presented for the fast and robust numerical solution of the 3D Poisson-Boltzmann Equation (PBE) governing the electrostatic interactions of large-scale biomolecules and highly charged multi-biomolecular assemblies such as ribosomes and viruses. The ACG offers numerous advantages over competing grid topologies such as regular 3D lattices and unstructured grids. For very large biological molecules and multi-biomolecule assemblies, the total number of grid-points is several orders of magnitude less than that required in a conventional lattice grid used in the current PBE solvers thus allowing the end user to obtain accurate and stable nonlinear PBE solutions on a desktop computer. Compared to tetrahedral-based unstructured grids, ACG offers a simpler hierarchical grid structure, which is naturally suited to multigrid, relieves indirect addressing requirements and uses fewer neighboring nodes in the finite difference stencils. Construction of the ACG and determination of the dielectric/ionic maps are straightforward, fast and require minimal user intervention. Charge singularities are eliminated by reformulating the problem to produce the reaction field potential in the molecular interior and the total electrostatic potential in the exterior ionic solvent region. This approach minimizes grid-dependency and alleviates the need for fine grid spacing near atomic charge sites. The technical portion of this paper contains three parts. First, the ACG and its construction for general biomolecular geometries are described. Next, a discrete approximation to the PBE upon this mesh is derived. Finally, the overall solution procedure and multigrid implementation are summarized. Results obtained with the ACG-based PBE solver are presented for: (i) a low dielectric spherical cavity, containing interior point charges, embedded in a high dielectric ionic solvent - analytical solutions are available for this case, thus allowing rigorous assessment of the solution accuracy; (ii) a pair of low dielectric charged spheres embedded in a ionic solvent to compute electrostatic interaction free energies as a function of the distance between sphere centers; (iii) surface potentials of proteins, nucleic acids and their larger-scale assemblies such as ribosomes; and (iv) electrostatic solvation free energies and their salt sensitivities - obtained with both linear and nonlinear Poisson-Boltzmann equation - for a large set of proteins. These latter results along with timings can serve as benchmarks for comparing the performance of different PBE solvers.

Molecular dynamics of ribosomal elongation factors G and Tu

Translation on the ribosome is controlled by external factors. During polypeptide lengthening, elongation factors EF-Tu and EF-G consecutively interact with the bacterial ribosome. EF-Tu binds and delivers an aminoacyl-tRNA to the ribosomal A site and EF-G helps translocate the tRNAs between their binding sites after the peptide bond is formed. These processes occur at the expense of GTP. EF-Tu:tRNA and EF-G are of similar shape, share a common binding site, and undergo large conformational changes on interaction with the ribosome. To characterize the internal motion of these two elongation factors, we used 25 ns long all-atom molecular dynamics simulations. We observed enhanced mobility of EF-G domains III, IV, and V and of tRNA in the EF-Tu:tRNA complex. EF-Tu:GDP complex acquired a configuration different from that found in the crystal structure of EF-Tu with a GTP analogue, showing conformational changes in the switch I and II regions. The calculated electrostatic properties of elongation factors showed no global similarity even though matching electrostatic surface patches were found around the domain I that contacts the ribosome, and in the GDP/GTP binding region.

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.

Non-bulk-like solvent behavior in the ribosome exit tunnel

As nascent proteins are synthesized by the ribosome, they depart via an exit tunnel running through the center of the large subunit. The exit tunnel likely plays an important part in various aspects of translation. Although water plays a key role in many bio-molecular processes, the nature of water confined to the exit tunnel has remained unknown. Furthermore, solvent in biological cavities has traditionally been characterized as either a continuous dielectric fluid, or a discrete tightly bound molecule. Using atomistic molecular dynamics simulations, we predict that the thermodynamic and kinetic properties of water confined within the ribosome exit tunnel are quite different from this simple two-state model. We find that the tunnel creates a complex microenvironment for the solvent resulting in perturbed rotational dynamics and heterogenous dielectric behavior. This gives rise to a very rugged solvation landscape and significantly retarded solvent diffusion. We discuss how this non-bulk-like solvent is likely to affect important biophysical processes such as sequence dependent stalling, co-translational folding, and antibiotic binding. We conclude with a discussion of the general applicability of these results to other biological cavities.

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.

Collective dynamics of the ribosomal tunnel revealed by elastic network modeling

The collective dynamics of the nascent polypeptide exit tunnel are investigated with the computationally efficient elastic network model using normal mode analysis. The calculated normal modes are considered individually and in linear combinations with different coefficients mimicking the phase angles between modes, in order to follow the mechanistic motions of tunnel wall residues. The low frequency fluctuations indicate three distinct regions along the tunnel-the entrance, the neck, and the exit-each having distinctly different domain motions. Generally, the lining of the entrance region moves in the exit direction, with the exit region having significantly larger motions, but in a perpendicular direction, whereas the confined neck region has rotational motions. Especially the universally conserved extensions of ribosomal proteins L4 and L22 located at the narrowest and mechanistically strategic region of tunnel undergo generally anti- or non-correlated motions, which may have an important role in nascent polypeptide gating mechanism. These motions appear to be sufficiently robust so as to be unaffected by the presence of a peptide chain in the tunnel.

Accurate prediction of protonation state as a prerequisite for reliable MM-PB(GB)SA binding free energy calculations of HIV-1 protease inhibitors

Binding free energies were calculated for the inhibitors lopinavir, ritonavir, saquinavir, indinavir, amprenavir, and nelfinavir bound to HIV-1 protease. An MMPB/SA-type analysis was applied to conformational samples from 3 ns explicit solvent molecular dynamics simulations of the enzyme-inhibitor complexes. Binding affinities and the sampled conformations of the inhibitor and enzyme were compared between different HIV-1 protease protonation states to find the most likely protonation state of the enzyme in the complex with each of the inhibitors. The resulting set of protonation states leads to good agreement between calculated and experimental binding affinities. Results from the MMPB/SA analysis are compared with an explicit/implicit hybrid scheme and with MMGB/SA methods. It is found that the inclusion of explicit water molecules may offer a slight advantage in reproducing absolute binding free energies while the use of the Generalized Born approximation significantly affects the accuracy of the calculated binding affinities.

Solvent reaction field potential inside an uncharged globular protein: a bridge between implicit and explicit solvent models?

The solvent reaction field potential of an uncharged protein immersed in simple point charge/extended explicit solvent was computed over a series of molecular dynamics trajectories, in total 1560 ns of simulation time. A finite, positive potential of 13-24 kbTec(-1) (where T=300 K), dependent on the geometry of the solvent-accessible surface, was observed inside the biomolecule. The primary contribution to this potential arose from a layer of positive charge density 1.0 A from the solute surface, on average 0.008 ec/A3, which we found to be the product of a highly ordered first solvation shell. Significant second solvation shell effects, including additional layers of charge density and a slight decrease in the short-range solvent-solvent interaction strength, were also observed. The impact of these findings on implicit solvent models was assessed by running similar explicit solvent simulations on the fully charged protein system. When the energy due to the solvent reaction field in the uncharged system is accounted for, correlation between per-atom electrostatic energies for the explicit solvent model and a simple implicit (Poisson) calculation is 0.97, and correlation between per-atom energies for the explicit solvent model and a previously published, optimized Poisson model is 0.99.

Correct and incorrect nucleotide incorporation pathways in DNA polymerase beta

Tracking the structural and energetic changes in the pathways of DNA replication and repair is central to the understanding of these important processes. Here we report favorable mechanisms of the polymerase-catalyzed phosphoryl transfer reactions corresponding to correct and incorrect nucleotide incorporations in the DNA by using a novel protocol involving energy minimizations, dynamics simulations, quasi-harmonic free energy calculations, and mixed quantum mechanics/molecular mechanics dynamics simulations. Though the pathway proposed may not be unique and invites variations, geometric and energetic arguments support the series of transient intermediates in the phosphoryl transfer pathways uncovered here for both the G:C and G:A systems involving a Grotthuss hopping mechanism of proton transfer between water molecules and the three conserved aspartate residues in pol beta's active-site. In the G:C system, the rate-limiting step is the initial proton hop with a free energy of activation of at least 17 kcal/mol, which corresponds closely to measured k(pol) values. Fidelity discrimination in pol beta can be explained by a significant loss of stability of the closed ternary complex of the enzyme in the G:A system and much higher activation energy of the initial step of nucleophilic attack, namely deprotonation of terminal DNA primer O3'H group. Thus, subtle differences in the enzyme active-site between matched and mismatched base pairs generate significant differences in catalytic performance.

Electrostatic properties of cowpea chlorotic mottle virus and cucumber mosaic virus capsids

Electrostatic properties of cowpea chlorotic mottle virus (CCMV) and cucumber mosaic virus (CMV) were investigated using numerical solutions to the Poisson-Boltzmann equation. Experimentally, it has been shown that CCMV particles swell in the absence of divalent cations when the pH is raised from 5 to 7. CMV, although structurally homologous, does not undergo this transition. An analysis of the calculated electrostatic potential confirms that a strong electrostatic repulsion at the calcium-binding sites in the CCMV capsid is most likely the driving force for the capsid swelling process during the release of calcium. The binding interaction between the encapsulated genome material (RNA) inside of the capsid and the inner capsid shell is weakened during the swelling transition. This probably aids in the RNA release process, but it is unlikely that the RNA is released through capsid openings due to unfavorable electrostatic interaction between the RNA and capsid inner shell residues at these openings. Calculations of the calcium binding energies show that Ca(2+) can bind both to the native and swollen forms of the CCMV virion. Favorable binding to the swollen form suggests that Ca(2+) ions can induce the capsid contraction and stabilize the native form.

SYMMETRY, FORM, AND SHAPE: Guiding Principles for Robustness in Macromolecular Machines

Computational studies of large macromolecular assemblages have come a long way during the past 10 years. With the explosion of computer power and parallel computing, timescales of molecular dynamics simulations have been extended far beyond the hundreds of picoseconds timescale. However, limitations remain for studies of large-scale conformational changes occurring on timescales beyond nanoseconds, especially for large macromolecules. In this review, we describe recent methods based on normal mode analysis that have enabled us to study dynamics on the microsecond timescale for large macromolecules using different levels of coarse graining, from atomically detailed models to those employing only low-resolution structural information. Emerging from such studies is a control principle for robustness in Nature's machines. We discuss this idea in the context of large-scale functional reorganization of the ribosome, virus particles, and the muscle protein myosin.

The electrostatic character of the ribosomal surface enables extraordinarily rapid target location by ribotoxins

Alpha-sarcin ribotoxins comprise a unique family of ribonucleases that cripple the ribosome by catalyzing endoribonucleolytic cleavage of ribosomal RNA at a specific location in the sarcin/ricin loop (SRL). The SRL structure alone is cleaved site-specifically by the ribotoxin, but the ribosomal context enhances the reaction rate by several orders of magnitude. We show that, for the alpha-sarcin-like ribotoxin restrictocin, this catalytic advantage arises from favorable electrostatic interactions with the ribosome. Restrictocin binds at many sites on the ribosomal surface and under certain conditions cleaves the SRL with a second-order rate constant of 1.7 x 10(10) M(-1) s(-1), a value that matches the predicted frequency of random restrictocin-ribosome encounters. The results suggest a mechanism of target location whereby restrictocin encounters ribosomes randomly and diffuses within the ribosomal electrostatic field to the SRL. These studies show a role for electrostatics in protein-ribosome recognition.

The role of positively charged amino acids and electrostatic interactions in the complex of U1A protein and U1 hairpin II RNA

Previous kinetic investigations of the N-terminal RNA recognition motif (RRM) domain of spliceosomal protein U1A, interacting with its RNA target U1 hairpin II, provided experimental evidence for a ‘lure and lock’ model of binding in which electrostatic interactions first guide the RNA to the protein, and close range interactions then lock the two molecules together. To further investigate the ‘lure’ step, here we examined the electrostatic roles of two sets of positively charged amino acids in U1A that do not make hydrogen bonds to the RNA: Lys20, Lys22 and Lys23 close to the RNA-binding site, and Arg7, Lys60 and Arg70, located on ‘top’ of the RRM domain, away from the RNA. Surface plasmon resonance-based kinetic studies, supplemented with salt dependence experiments and molecular dynamics simulation, indicate that Lys20 predominantly plays a role in association, while nearby residues Lys22 and Lys23 appear to be at least as important for complex stability. In contrast, kinetic analyses of residues away from the RNA indicate that they have a minimal effect on association and stability. Thus, well-positioned positively charged residues can be important for both initial complex formation and complex maintenance, illustrating the multiple roles of electrostatic interactions in protein–RNA complexes.

The Extent of Cooperativity of Protein Motions Observed with Elastic Network Models Is Similar for Atomic and Coarser-Grained Models

Coarse-grained elastic network models have been successful in determining functionally relevant collective motions. The level of coarse-graining, however, has usually focused on the level of one point per residue. In this work, we compare the applicability of elastic network models over a broader range of representational scales. We apply normal mode analysis for multiple scales on a high-resolution protein data set using various cutoff radii to define the residues considered to be interacting, or the extent of cooperativity of their motions. These scales include the residue-, atomic-, proton-, and explicit solvent-levels. Interestingly, atomic, proton, and explicit solvent level calculations all provide similar results at the same cutoff value, with the computed mean-square fluctuations showing only a slightly higher correlation (0.61) with the experimental temperature factors from crystallography than the results of the residue-level coarse-graining. The qualitative behavior of each level of coarse graining is similar at different cutoff values. The correlations between these fluctuations and the number of internal contacts improve with increased cutoff values. Our results demonstrate that atomic level elastic network models provide an improved representation for the collective motions of proteins compared to the coarse-grained models.

Simulating movement of tRNA into the ribosome during decoding

Decoding is the key step during protein synthesis that enables information transfer from RNA to protein, a process critical for the survival of all organisms. We have used large-scale (2.64 x 10(6) atoms) all-atom simulations of the entire ribosome to understand a critical step of decoding. Although the decoding problem has been studied for more than four decades, the rate-limiting step of cognate tRNA selection has only recently been identified. This step, known as accommodation, involves the movement inside the ribosome of the aminoacyl-tRNA from the partially bound "A/T" state to the fully bound "A/A" state. Here, we show that a corridor of 20 universally conserved ribosomal RNA bases interacts with the tRNA during the accommodation movement. Surprisingly, the tRNA is impeded by the A-loop (23S helix 92), instead of enjoying a smooth transition to the A/A state. In particular, universally conserved 23S ribosomal RNA bases U2492, C2556, and C2573 act as a 3D gate, causing the acceptor stem to pause before allowing entrance into the peptidyl transferase center. Our simulations demonstrate that the flexibility of the acceptor stem of the tRNA, in addition to flexibility of the anticodon arm, is essential for tRNA selection. This study serves as a template for simulating conformational changes in large (>10(6) atoms) biological and artificial molecular machines.

Exploring global motions and correlations in the ribosome

We studied slower global coupled motions of the ribosome with half a microsecond of coarse-grained molecular dynamics. A low-resolution anharmonic network model that allows for the evolution of tertiary structure and long-scale sampling was developed and parameterized. Most importantly, we find that functionally important movements of L7/L12 and L1 lateral stalks are anticorrelated. Other principal directions of motions include widening of the tRNA cleft and the rotation of the small subunit which occurs as one block and is in phase with the movement of L1 stalk. The effect of the dynamical correlation pattern on the elongation process is discussed. Small fluctuations of the 3' tRNA termini and anticodon nucleotides show tight alignment of substrates for the reaction. Our model provides an efficient and reliable way to study the dynamics of large biomolecular systems composed of both proteins and nucleic acids.

Electrostatic interaction between RNA and protein capsid in cowpea chlorotic mottle virus simulated by a coarse-grain RNA model and a Monte Carlo approach

Although many viruses have been crystallized and the protein capsid structures have been determined by x-ray crystallography, the nucleic acids often cannot be resolved. This is especially true for RNA viruses. The lack of information about the conformation of DNA/RNA greatly hinders our understanding of the assembly mechanism of various viruses. Here we combine a coarse-grain model and a Monte Carlo method to simulate the distribution of viral RNA inside the capsid of cowpea chlorotic mottle virus. Our results show that there is very strong interaction between the N-terminal residues of the capsid proteins, which are highly positive charged, and the viral RNA. Without these residues, the binding energy disfavors the binding of RNA by the capsid. The RNA forms a shell close to the capsid with the highest densities associated with the capsid dimers. These high-density regions are connected to each other in the shape of a continuous net of triangles. The overall icosahedral shape of the net overlaps with the capsid subunit icosahedral organization. Medium density of RNA is found under the pentamers of the capsid. These findings are consistent with experimental observations.