Recognition of the regulatory nascent chain TnaC by the ribosome

A machine learning approach uncovers principles and determinants of eukaryotic ribosome pausing

Nonuniform local translation speed dictates diverse protein biogenesis outcomes. To unify known and uncover unknown principles governing eukaryotic elongation rate, we developed a machine learning pipeline to analyze RiboSeq datasets. We find that the chemical nature of the incoming amino acid determines how codon optimality influences elongation rate, with hydrophobic residues more dependent on transfer RNA (tRNA) levels than charged residues. Unexpectedly, we find that wobble interactions exert a widespread effect on elongation pausing, with wobble-mediated decoding being slower than Watson-Crick decoding, irrespective of tRNA levels. Applying our ribosome pausing principles to ribosome collisions reveals that disomes arise upon apposition of fast-decoding and slow-decoding signatures. We conclude that codon choice and tRNA pools are evolutionarily constrained to harmonize elongation rate with cotranslational folding while minimizing wobble pairing and deleterious stalling.

The energy landscape of the ribosome

The ribosome is a prototypical assembly that can be used to establish general principles and techniques for the study of biological molecular machines. Motivated by the fact that the dynamics of every biomolecule is governed by an underlying energy landscape, there has been great interest to understand and quantify ribosome energetics. In the present review, we will focus on theoretical and computational strategies for probing the interactions that shape the energy landscape of the ribosome, with an emphasis on more recent studies of the elongation cycle. These efforts include the application of quantum mechanical methods for describing chemical kinetics, as well as classical descriptions to characterize slower (microsecond to millisecond) large-scale (10-100 Å) rearrangements, where motion is described in terms of diffusion across an energy landscape. Together, these studies provide broad insights into the factors that control a diverse range of dynamical processes in this assembly.

Atomistic simulations of the Escherichia coli ribosome provide selection criteria for translationally active substrates

As genetic code expansion advances beyond L-α-amino acids to backbone modifications and new polymerization chemistries, delineating what substrates the ribosome can accommodate remains a challenge. The Escherichia coli ribosome tolerates non-L-α-amino acids in vitro, but few structural insights that explain how are available, and the boundary conditions for efficient bond formation are so far unknown. Here we determine a high-resolution cryogenic electron microscopy structure of the E. coli ribosome containing α-amino acid monomers and use metadynamics simulations to define energy surface minima and understand incorporation efficiencies. Reactive monomers across diverse structural classes favour a conformational space where the aminoacyl-tRNA nucleophile is <4 Å from the peptidyl-tRNA carbonyl with a Bürgi-Dunitz angle of 76-115°. Monomers with free energy minima that fall outside this conformational space do not react efficiently. This insight should accelerate the in vivo and in vitro ribosomal synthesis of sequence-defined, non-peptide heterooligomers.

Engineering Toehold-Mediated Switches for Native RNA Detection and Regulation in Bacteria

RNA switches are versatile tools in synthetic biology for sensing and regulation applications. The discoveries of RNA-mediated translational and transcriptional control have facilitated the development of complexde novodesigns of RNA switches. Specifically, RNA toehold-mediated switches, in which binding to the toehold sensing domain controls the transition between switch states via strand displacement, have been extensively adapted for coupling systems responses to specifictrans-RNA inputs. This review highlights some of the challenges associated with applying these switches for native RNA detectionin vivo, including transferability between organisms. The applicability and design considerations of toehold-mediated switches are discussed by highlighting twelve recently developed switch designs. This review finishes with future perspectives to address current gaps in the field, particularly regarding the power of structural prediction algorithms for improved in vivo functionality of RNA switches.

Exploring Allosteric Signaling in the Exit Tunnel of the Bacterial Ribosome by Molecular Dynamics Simulations and Residue Network Model

The bacterial ribosomal tunnel is equipped with numerous sites highly sensitive to the course of the translation process. This study investigates allosteric pathways linking distant functional sites that collaboratively play a role either in translation regulation or recruitment of chaperones. We apply perturbation response scanning (PRS) analysis to 700 ns long and 500 ns long coarse-grained molecular dynamics simulations of E. coli and T. thermophilus large subunits, respectively, to reveal nucleotides/residues with the ability to transmit perturbations by dynamic rationale. We also use the residue network model with the k-shortest pathways method to calculate suboptimal pathways based on the contact topology of the ribosomal tunnel of E. coli crystal structure and 101 ClustENM generated conformers of T. thermophilus large subunit. In the upper part of the tunnel, results suggest that A2062 and A2451 can communicate in both directions for translation stalling, mostly through dynamically coupled C2063, C2064, and A2450. For a similar purpose, U2585 and U2586 are coupled with A2062, while they are also sensitive to uL4 and uL22 at the constriction region through two different pathways at the opposite sides of the tunnel wall. In addition, the constriction region communicates with the chaperone binding site on uL23 at the solvent side but through few nucleotides. Potential allosteric communication pathways between the lower part of the tunnel and chaperone binding site mostly use the flexible loop of uL23, while A1336–G1339 provide a suboptimal pathway. Both species seem to employ similar mechanisms in the long tunnel, where a non-conserved cavity at the bacterial uL23 and 23S rRNA interface is proposed as a novel drug target.

Large-scale simulations of nucleoprotein complexes: ribosomes, nucleosomes, chromatin, chromosomes and CRISPR

Recent advances in biotechnology such as Hi-C, CRISPR/Cas9 and ribosome display have placed nucleoprotein complexes at center stage. Understanding the structural dynamics of these complexes aids in optimizing protocols and interpreting data for these new technologies. The integration of simulation and experiment has helped advance mechanistic understanding of these systems. Coarse-grained simulations, reduced-description models, and explicit solvent molecular dynamics simulations yield useful complementary perspectives on nucleoprotein complex structural dynamics. When combined with Hi-C, cryo-EM, and single molecule measurements, these simulations integrate disparate forms of experimental data into a coherent mechanism.

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

With both catalytic and genetic functions, ribonucleic acid (RNA) is perhaps the most pluripotent chemical species in molecular biology, and its functions are intimately linked to its structure and dynamics. Computer simulations, and in particular atomistic molecular dynamics (MD), allow structural dynamics of biomolecular systems to be investigated with unprecedented temporal and spatial resolution. We here provide a comprehensive overview of the fast-developing field of MD simulations of RNA molecules. We begin with an in-depth, evaluatory coverage of the most fundamental methodological challenges that set the basis for the future development of the field, in particular, the current developments and inherent physical limitations of the atomistic force fields and the recent advances in a broad spectrum of enhanced sampling methods. We also survey the closely related field of coarse-grained modeling of RNA systems. After dealing with the methodological aspects, we provide an exhaustive overview of the available RNA simulation literature, ranging from studies of the smallest RNA oligonucleotides to investigations of the entire ribosome. Our review encompasses tetranucleotides, tetraloops, a number of small RNA motifs, A-helix RNA, kissing-loop complexes, the TAR RNA element, the decoding center and other important regions of the ribosome, as well as assorted others systems. Extended sections are devoted to RNA-ion interactions, ribozymes, riboswitches, and protein/RNA complexes. Our overview is written for as broad of an audience as possible, aiming to provide a much-needed interdisciplinary bridge between computation and experiment, together with a perspective on the future of the field.

The role of recombination in evolutionary adaptation of Escherichia coli to a novel nutrient

The benefits and detriments of recombination for adaptive evolution have been studied both theoretically and experimentally, with conflicting predictions and observations. Most pertinent experiments examine recombination's effects in an unchanging environment and do not study its genomewide effects. Here, we evolved six replicate populations of either highly recombining R⁺ or lowly recombining R^- E. coli strains in a changing environment, by introducing the novel nutrients L-arabinose or indole into the environment. The experiment's ancestral strains are not viable on these nutrients, but 130 generations of adaptive evolution were sufficient to render them viable. Recombination conferred a more pronounced advantage to populations adapting to indole. To study the genomic changes associated with this advantage, we sequenced the genomes of 384 clones isolated from selected replicates at the end of the experiment. These genomes harbour complex changes that range from point mutations to large-scale DNA amplifications. Among several candidate adaptive mutations, those in the tryptophanase regulator tnaC stand out, because the tna operon in which it resides has a known role in indole metabolism. One of the highly recombining populations also shows a significant excess of large-scale segmental DNA amplifications that include the tna operon. This lineage also shows a unique and potentially adaptive combination of point mutations and DNA amplifications that may have originated independently from one another, to be joined later by recombination. Our data illustrate that the advantages of recombination for adaptive evolution strongly depend on the environment and that they can be associated with complex genomic changes.

Effect of Nascent Peptide Steric Bulk on Elongation Kinetics in the Ribosome Exit Tunnel

All proteins are synthesized by the ribosome, a macromolecular complex that accomplishes the life-sustaining tasks of faithfully decoding mRNA and catalyzing peptide bond formation at the peptidyl transferase center (PTC). The ribosome has evolved an exit tunnel to host the elongating new peptide, protect it from proteolytic digestion, and guide its emergence. It is here that the nascent chain begins to fold. This folding process depends on the rate of translation at the PTC. We report here that besides PTC events, translation kinetics depend on steric constraints on nascent peptide side chains and that confined movements of cramped side chains within and through the tunnel fine-tune elongation rates.

The molecular choreography of protein synthesis: translational control, regulation, and pathways

Translation of proteins by the ribosome regulates gene expression, with recent results underscoring the importance of translational control. Misregulation of translation underlies many diseases, including cancer and many genetic diseases. Decades of biochemical and structural studies have delineated many of the mechanistic details in prokaryotic translation, and sketched the outlines of eukaryotic translation. However, translation may not proceed linearly through a single mechanistic pathway, but likely involves multiple pathways and branchpoints. The stochastic nature of biological processes would allow different pathways to occur during translation that are biased by the interaction of the ribosome with other translation factors, with many of the steps kinetically controlled. These multiple pathways and branchpoints are potential regulatory nexus, allowing gene expression to be tuned at the translational level. As research focus shifts toward eukaryotic translation, certain themes will be echoed from studies on prokaryotic translation. This review provides a general overview of the dynamic data related to prokaryotic and eukaryotic translation, in particular recent findings with single-molecule methods, complemented by biochemical, kinetic, and structural findings. We will underscore the importance of viewing the process through the viewpoints of regulation, translational control, and heterogeneous pathways.

Regional discrimination and propagation of local rearrangements along the ribosomal exit tunnel

All proteins, from bacteria to man, are made in the ribosome and are elongated, one residue at a time, at the peptidyl transferase center. This growing peptide chain wends its way through the ribosomal tunnel to the exit port, ~100Å from the peptidyl transferase center. We have identified locations in the tunnel that sense and respond to single side chains of the nascent peptide to induce local conformational changes. Moreover, side-chain sterics and rearrangements deep in the tunnel influence the disposition of residues 45Å away at the exit port and are consistent with side-chain-induced axial retraction of the peptide backbone. These coupled responses are neither haphazard nor uniform along the tunnel. Rather, they are confined to discriminating zones in the tunnel and are sequence specific. Such discerning communication may contribute to folding events and mechanisms governing sequence-specific signaling between different regions of the tunnel during translation.

Conjugates of amino acids and peptides with 5-o-mycaminosyltylonolide and their interaction with the ribosomal exit tunnel

During protein synthesis the nascent polypeptide chain (NC) extends through the ribosomal exit tunnel (NPET). Also, the large group of macrolide antibiotics binds in the nascent peptide exit tunnel. In some cases interaction of NC with NPET leads to the ribosome stalling, a significant event in regulation of translation. In other cases NC-ribosome interactions lead to pauses in translation that play an important role in cotranslational folding of polypeptides emerging from the ribosome. The precise mechanism of NC recognition in NPET as well as factors that determine NC conformation in the ribosomal tunnel are unknown. A number of derivatives of the macrolide antibiotic 5-O-mycaminosyltylonolide (OMT) containing N-acylated amino acid or peptide residues were synthesized in order to study potential sites of NC-NPET interactions. The target compounds were prepared by conjugation of protected amino acids and peptides with the C23 hydroxyl group of the macrolide. These OMT derivatives showed high although varying abilities to inhibit the firefly luciferase synthesis in vitro. Three glycil-containing derivatives appeared to be strong inhibitors of translation, more potent than parental OMT. Molecular dynamics (MD) simulation of complexes of tylosin, OMT, and some of OMT derivatives with the large ribosomal subunit of E. coli illuminated a plausible reason for the high inhibitory activity of Boc-Gly-OMT. In addition, the MD study detected a new putative site of interaction of the nascent polypeptide chain with the NPET walls.

Interactions of the TnaC nascent peptide with rRNA in the exit tunnel enable the ribosome to respond to free tryptophan

A transcriptional attenuation mechanism regulates expression of the bacterial tnaCAB operon. This mechanism requires ribosomal arrest induced by the regulatory nascent TnaC peptide in response to free L-tryptophan (L-Trp). In this study we demonstrate, using genetic and biochemical analyses, that in Escherichia coli, TnaC residue I19 and 23S rRNA nucleotide A2058 are essential for the ribosome's ability to sense free L-Trp. We show that the mutational change A2058U in 23S rRNA reduces the concentration dependence of L-Trp-mediated tna operon induction, whereas the TnaC I19L change suppresses this phenotype, restoring the sensitivity of the translating A2058U mutant ribosome to free L-Trp. These findings suggest that interactions between TnaC residue I19 and 23S rRNA nucleotide A2058 contribute to the creation of a regulatory L-Trp binding site within the translating ribosome.

Single-molecule nanometry for biological physics

Precision measurement is a hallmark of physics but the small length scale (∼nanometer) of elementary biological components and thermal fluctuations surrounding them challenge our ability to visualize their action. Here, we highlight the recent developments in single-molecule nanometry where the position of a single fluorescent molecule can be determined with nanometer precision, reaching the limit imposed by the shot noise, and the relative motion between two molecules can be determined with ∼0.3 nm precision at ∼1 ms time resolution, as well as how these new tools are providing fundamental insights into how motor proteins move on cellular highways. We will also discuss how interactions between three and four fluorescent molecules can be used to measure three and six coordinates, respectively, allowing us to correlate the movements of multiple components. Finally, we will discuss recent progress in combining angstrom-precision optical tweezers with single-molecule fluorescent detection, opening new windows for multi-dimensional single-molecule nanometry for biological physics.

Cryo-electron microscopy modeling by the molecular dynamics flexible fitting method

The increasing power and popularity of cryo-electron microscopy (cryo-EM) in structural biology brought about the development of so-called hybrid methods, which permit the interpretation of cryo-EM density maps beyond their nominal resolution in terms of atomic models. The Cryo-EM Modeling Challenge 2010 is the first community effort to bring together developers of hybrid methods as well as cryo-EM experimentalists. Participating in the challenge, the molecular dynamics flexible fitting (MDFF) method was applied to a number of cryo-EM density maps. The results are described here with special emphasis on the use of symmetry-based restraints to improve the quality of atomic models derived from density maps of symmetric complexes; on a comparison of the stereochemical quality of atomic models resulting from different hybrid methods; and on application of MDFF to electron crystallography data.

Mechanisms of SecM-mediated stalling in the ribosome

Nascent-peptide modulation of translation is a common regulatory mechanism of gene expression. In this mechanism, while the nascent peptide is still in the exit tunnel of the ribosome, it induces translational pausing, thereby controlling the expression of downstream genes. One example is SecM, which inhibits peptide-bond formation in the ribosome's peptidyl transferase center (PTC) during its own translation, upregulating the expression of the protein translocase SecA. Although biochemical experiments and cryo-electron microscopy data have led to the identification of some residues involved in SecM recognition, the full pathway of interacting residues that connect SecM to the PTC through the ribosome has not yet been conclusively established. Here, using the cryo-electron microscopy data, we derived the first (to our knowledge) atomic model of the SecM-stalled ribosome via molecular-dynamics flexible fitting, complete with P- and A-site tRNAs. Subsequently, we carried out simulations of native and mutated SecM-stalled ribosomes to investigate possible interaction pathways between a critical SecM residue, R163, and the PTC. In particular, the simulations reveal the role of SecM in altering the position of the tRNAs in the ribosome, and thus demonstrate how the presence of SecM in the exit tunnel induces stalling. Finally, steered molecular-dynamics simulations in which SecM was pulled toward the tunnel exit suggest how SecA interacting with SecM from outside the ribosome relieves stalling.

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.

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.

Arginine changes the conformation of the arginine attenuator peptide relative to the ribosome tunnel

The fungal arginine attenuator peptide (AAP) is a regulatory peptide that controls ribosome function. As a nascent peptide within the ribosome exit tunnel, it acts to stall ribosomes in response to arginine (Arg). We used three approaches to probe the molecular basis for stalling. First, PEGylation assays revealed that the AAP did not undergo overall compaction in the tunnel in response to Arg. Second, site-specific photocross-linking showed that Arg altered the conformation of the wild-type AAP, but not of nonfunctional mutants, with respect to the tunnel. Third, using time-resolved spectral measurements with a fluorescent probe placed in the nascent AAP, we detected sequence-specific changes in the disposition of the AAP near the peptidyltransferase center in response to Arg. These data provide evidence that an Arg-induced change in AAP conformation and/or environment in the ribosome tunnel is important for stalling.

Crucial elements that maintain the interactions between the regulatory TnaC peptide and the ribosome exit tunnel responsible for Trp inhibition of ribosome function

Translation of the TnaC nascent peptide inhibits ribosomal activity in the presence of l-tryptophan, inducing expression of the tnaCAB operon in Escherichia coli. Using chemical methylation, this work reveals how interactions between TnaC and the ribosome are affected by mutations in both molecules. The presence of the TnaC-tRNA^Pro peptidyl-tRNA within the ribosome protects the 23S rRNA nucleotide U2609 against chemical methylation. Such protection was not observed in mutant ribosomes containing changes in 23S rRNA nucleotides of the A748–A752 region. Nucleotides A752 and U2609 establish a base-pair interaction. Most replacements of either A752 or U2609 affected Trp induction of a TnaC-regulated LacZ reporter. However, the single change A752G, or the dual replacements A752G and U2609C, maintained Trp induction. Replacements at the conserved TnaC residues W12 and D16 also abolished the protection of U2609 by TnaC-tRNA^Pro against chemical methylation. These data indicate that the TnaC nascent peptide in the ribosome exit tunnel interacts with the U2609 nucleotide when the ribosome is Trp responsive. This interaction is affected by mutational changes in exit tunnel nucleotides of 23S rRNA, as well as in conserved TnaC residues, suggesting that they affect the structure of the exit tunnel and/or the nascent peptide configuration in the tunnel.

Ribosome assembly factors prevent premature translation initiation by 40S assembly intermediates

Ribosome assembly in eukaryotes requires approximately 200 essential assembly factors (AFs) and occurs through ordered events that initiate in the nucleolus and culminate in the cytoplasm. Here, we present the electron cryo-microscopy (cryo-EM) structure of a late cytoplasmic 40S ribosome assembly intermediate from Saccharomyces cerevisiae at 18 angstrom resolution. We obtained cryo-EM reconstructions of preribosomal complexes lacking individual components to define the positions of all seven AFs bound to this intermediate. These late-binding AFs are positioned to prevent each step in the translation initiation pathway. Together, they obstruct the binding sites for initiation factors, prevent the opening of the messenger RNA channel, block 60S subunit joining, and disrupt the decoding site. These redundant mechanisms probably ensure that pre-40S particles do not enter the translation pathway, which would result in their rapid degradation.

Nascent peptide side chains induce rearrangements in distinct locations of the ribosomal tunnel

Although we have numerous structures of ribosomes, none disclose side-chain rearrangements of the nascent peptide during chain elongation. This study reports for the first time that rearrangement of the peptide and/or tunnel occurs in distinct regions of the tunnel and is directed by the unique primary sequence of each nascent peptide. In the tunnel mid-region, the accessibility of an introduced cysteine to a series of novel hydrophilic maleimide reagents increases with increasing volume of the adjacent chain residue, a sensitivity not manifest at the constriction and exit port. This surprising result reveals molecular movements not yet resolvable from structural studies. These findings map solvent-accessible volumes along the tunnel and provide novel insights critical to our understanding of allosteric communication within the ribosomal tunnel, translational arrest, chaperone interaction, folding, and rates of elongation.

Free energy of nascent-chain folding in the translocon

During their synthesis, many water-soluble proteins and nearly all membrane proteins transit through a protein-conducting channel in the membrane, the Sec translocon, from where they are inserted into the lipid bilayer. Increasing evidence indicates that folding of the nascent protein begins already within the ribosomal exit tunnel in a sequence- and environment-dependent fashion. To examine the effects of the translocon on the nascent-chain folding, we have calculated the potential of mean force for α-helix formation of a 10-alanine oligopeptide as a function of its position within the translocon channel. We find that the predominant conformational states, α-helical and extended, reflect those found for the peptide in water. However, the translocon, via its surface properties and its variable diameter, shifts the equilibrium in favor of the α-helical state. Thus, we suggest that the translocon facilitates not only the insertion of membrane proteins into the bilayer but also their folding.

Nascent polypeptide sequences that influence ribosome function

Ribosomes catalyze protein synthesis using transfer RNAs and auxiliary proteins. Historically, ribosomes have been considered nonspecific translational machines, having no regulatory functions. However, a new class of regulatory mechanisms has been discovered that is based on interactions occurring within the ribosomal peptide exit tunnel that result in ribosome stalling during translation of an appropriate mRNA segment. These discoveries reveal an unexpectedly dynamic role ribosomes play in regulating their own activity. By using nascent leader peptides in combination with bound specific amino acids or antibiotics, ribosome functions can be altered significantly resulting in regulated expression of downstream coding regions. This review summarizes relevant findings in recent articles and outlines our current understanding of nascent peptide-induced ribosome stalling in regulating gene expression.

Connecting energy landscapes with experimental rates for aminoacyl-tRNA accommodation in the ribosome

Using explicit-solvent simulations of the 70S ribosome, the barrier-crossing attempt frequency was calculated for aminoacyl-tRNA elbow-accommodation. In seven individual trajectories (200-300 ns, each, for an aggregate time of 2.1 μs), the relaxation time of tRNA structural fluctuations was determined to be ∼10 ns, and the barrier-crossing attempt frequency of tRNA accommodation is ∼1-10 μs(-1). These calculations provide a quantitative relationship between the free-energy barrier and experimentally measured rates of accommodation, which demonstrate that the free-energy barrier of elbow-accommodation is less than 15 k(B)T, in vitro and in vivo.

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.

Applications of the molecular dynamics flexible fitting method

In recent years, cryo-electron microscopy (cryo-EM) has established itself as a key method in structural biology, permitting the structural characterization of large biomolecular complexes in various functional states. The data obtained through single-particle cryo-EM has recently seen a leap in resolution thanks to landmark advances in experimental and computational techniques, resulting in sub-nanometer resolution structures being obtained routinely. The remaining gap between these data and revealing the mechanisms of molecular function can be closed through hybrid modeling tools that incorporate known atomic structures into the cryo-EM data. One such tool, molecular dynamics flexible fitting (MDFF), uses molecular dynamics simulations to combine structures from X-ray crystallography with cryo-EM density maps to derive atomic models of large biomolecular complexes. The structures furnished by MDFF can be used subsequently in computational investigations aimed at revealing the dynamics of the complexes under study. In the present work, recent applications of MDFF are presented, including the interpretation of cryo-EM data of the ribosome at different stages of translation and the structure of a membrane-curvature-inducing photosynthetic complex.

Ionic partition and transport in multi-ionic channels: a molecular dynamics simulation study of the OmpF bacterial porin

We performed all-atom molecular dynamics simulations studying the partition of ions and the ionic current through the bacterial porin OmpF and two selected mutants. The study is motivated by new, interesting experimental findings concerning their selectivity and conductance behavior at neutral pH. The mutations considered here are designed to study the effect of removal of negative charges present in the constriction zone of the wild-type OmpF channel (which contains, on one side, a cluster with three positive residues, and on the other side, two negatively charged residues). Our results show that these mutations induce an exclusion of cations from the constriction zone of the channel, substantially reducing the flow of cations. In fact, the partition of ions inside the mutant channels is strongly inhomogeneous, with regions containing an excess of cations and regions containing an excess of anions. Interestingly, the overall number of cations inside the channel is larger than the number of anions, this excess being different for each protein channel. We found that the differences in ionic charge inside these channels are justified by the differences in electric charge between the wild-type OmpF and the mutants, following an electroneutral balance.

GPU-accelerated molecular modeling coming of age

Graphics processing units (GPUs) have traditionally been used in molecular modeling solely for visualization of molecular structures and animation of trajectories resulting from molecular dynamics simulations. Modern GPUs have evolved into fully programmable, massively parallel co-processors that can now be exploited to accelerate many scientific computations, typically providing about one order of magnitude speedup over CPU code and in special cases providing speedups of two orders of magnitude. This paper surveys the development of molecular modeling algorithms that leverage GPU computing, the advances already made and remaining issues to be resolved, and the continuing evolution of GPU technology that promises to become even more useful to molecular modeling. Hardware acceleration with commodity GPUs is expected to benefit the overall computational biology community by bringing teraflops performance to desktop workstations and in some cases potentially changing what were formerly batch-mode computational jobs into interactive tasks.