Statics of the ribosomal exit tunnel: implications for cotranslational peptide folding, elongation regulation, and antibiotics binding

Design, synthesis and biological evaluation of amphiphilic benzopyran derivatives as potent antibacterial agents against multidrug-resistant bacteria

Antimicrobial resistance has emerged as a significant threat to global public health. To develop novel, high efficiency antibacterial alternatives to combat multidrug-resistant bacteria, A total of thirty-two novel amphiphilic benzopyran derivatives by mimicking the structure and function of antimicrobial peptides were designed and synthesized. Among them, the most promising compounds 4h and 17e displayed excellent antibacterial activity against Gram-positive bacteria (MICs = 1-4 μg/mL) with weak hemolytic activity and good membrane selectivity. Additionally, compounds 4h and 17e had rapid bactericidal properties, low resistance frequency, good plasma stability, and strong capabilities of inhibiting and eliminating bacterial biofilms. Mechanistic studies revealed that compounds 4h and 17e could effectively disrupt the integrity of bacterial cell membranes, and accompanied by an increase in intracellular reactive oxygen species and the leakage of proteins and DNA, ultimately leading to bacterial death. Notably, compound 4h exhibited comparable in vivo antibacterial potency in a mouse septicemia model infected by Staphylococcus aureus ATCC43300, as compared to vancomycin. These findings indicated that 4h might be a promising antibacterial candidate to combat antimicrobial resistance.

Macrolide therapy in Pseudomonas aeruginosa infections causes uL4 ribosomal protein mutations leading to high-level resistance

OBJECTIVES

Pseudomonas aeruginosa colonizes the cystic fibrosis (CF) airways causing chronic bacterial lung infections. CF patients are routinely treated with macrolides, however, P. aeruginosa is considered insusceptible as consequence of inadequate susceptibility testing leaving resistance mechanism completely overlooked. Here, we investigated a new mechanism of macrolide resistance caused by ribosomal protein mutations.

METHODS

Investigating a longitudinal collection of 529 isolates from CF patients and analysing 5758 protein sequences from different sources, mutations in P. aeruginosa's ribosomal proteins connected to macrolide resistance were identified. Using a modified susceptibility testing protocol, isolates harbouring a mutated uL4 ribosomal protein were tested for resistance against macrolide antibiotics and macrolide-induced quorum sensing modulation. Proteome and ribosome profiling were applied to assess the impact of the mutations on the bacterial physiology.

RESULTS

Five uL4 mutations were identified in isolates from different CF patients. Most mapped to the conserved loop region of uL4 and resulted in increased macrolide tolerance (>10-fold relative to wt strains). Greater concentrations (>10-fold) of macrolide antibiotic were needed to inhibit the growth, reduce swimming motility, and induce redox sensitivity of the uL4 mutants. 16 proteins involved in ribosome adaptation displayed altered expression possibly to compensate for the uL4 mutations, which changed the ribosome stoichiometry without negatively affecting bacterial physiology.

CONCLUSIONS

Macrolide antibiotics should, therefore, be considered as active antimicrobial agents against P. aeruginosa and resistance development should be contemplated when patients are treated with prolonged courses of macrolides. Importantly, improved macrolide susceptibility testing is necessary for the detection of resistant bacteria.

An intrinsically disordered nascent protein interacts with specific regions of the ribosomal surface near the exit tunnel

The influence of the ribosome on nascent chains is poorly understood, especially in the case of proteins devoid of signal or arrest sequences. Here, we provide explicit evidence for the interaction of specific ribosomal proteins with ribosome-bound nascent chains (RNCs). We target RNCs pertaining to the intrinsically disordered protein PIR and a number of mutants bearing a variable net charge. All the constructs analyzed in this work lack N-terminal signal sequences. By a combination chemical crosslinking and Western-blotting, we find that all RNCs interact with ribosomal protein L23 and that longer nascent chains also weakly interact with L29. The interacting proteins are spatially clustered on a specific region of the large ribosomal subunit, close to the exit tunnel. Based on chain-length-dependence and mutational studies, we find that the interactions with L23 persist despite drastic variations in RNC sequence. Importantly, we also find that the interactions are highly Mg+2-concentration-dependent. This work is significant because it unravels a novel role of the ribosome, which is shown to engage with the nascent protein chain even in the absence of signal or arrest sequences.

Recent Advances in the Development of Macrolide Antibiotics as Antimicrobial Agents

The chemical modification of natural products has been a major method in the discovery and synthesis of new macrolide antibiotics (MA) to treat a variety of infectious diseases. However, a lot of MA obtained in the above methods are no longer effective, because the bacteria quickly develop their resistance to these new macrolides, which has become a great threat to successful treatment of infectious diseases, such as infections of the respiratory system and urinary system. In this paper, total synthetic methods for MA that include erythromycin A (ERY), azithromycin (AZM), the clinical candidate solithromycin (CEM-101), as well as 14-membered and 15-membered azaketolides have been systematically reviewed on the basis of the literature reported previously. The total synthetic methods we describe here helps to accelerate the discovery of newer MA to deal with the serious problem of bacterial resistance.

Root of the Tree: The Significance, Evolution, and Origins of the Ribosome

The ribosome is an ancient molecular fossil that provides a telescope to the origins of life. Made from RNA and protein, the ribosome translates mRNA to coded protein in all living systems. Universality, economy, centrality and antiquity are ingrained in translation. The translation machinery dominates the set of genes that are shared as orthologues across the tree of life. The lineage of the translation system defines the universal tree of life. The function of a ribosome is to build ribosomes; to accomplish this task, ribosomes make ribosomal proteins, polymerases, enzymes, and signaling proteins. Every coded protein ever produced by life on Earth has passed through the exit tunnel, which is the birth canal of biology. During the root phase of the tree of life, before the last common ancestor of life (LUCA), exit tunnel evolution is dominant and unremitting. Protein folding coevolved with evolution of the exit tunnel. The ribosome shows that protein folding initiated with intrinsic disorder, supported through a short, primitive exit tunnel. Folding progressed to thermodynamically stable β-structures and then to kinetically trapped α-structures. The latter were enabled by a long, mature exit tunnel that partially offset the general thermodynamic tendency of all polypeptides to form β-sheets. RNA chaperoned the evolution of protein folding from the very beginning. The universal common core of the ribosome, with a mass of nearly 2 million Daltons, was finalized by LUCA. The ribosome entered stasis after LUCA and remained in that state for billions of years. Bacterial ribosomes never left stasis. Archaeal ribosomes have remained near stasis, except for the superphylum Asgard, which has accreted rRNA post LUCA. Eukaryotic ribosomes in some lineages appear to be logarithmically accreting rRNA over the last billion years. Ribosomal expansion in Asgard and Eukarya has been incremental and iterative, without substantial remodeling of pre-existing basal structures. The ribosome preserves information on its history.

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.

Structural insights of lincosamides targeting the ribosome of Staphylococcus aureus

Antimicrobial resistance within a wide range of pathogenic bacteria is an increasingly serious threat to global public health. Among these pathogenic bacteria are the highly resistant, versatile and possibly aggressive bacteria, Staphylococcus aureus. Lincosamide antibiotics were proved to be effective against this pathogen. This small, albeit important group of antibiotics is mostly active against Gram-positive bacteria, but also used against selected Gram-negative anaerobes and protozoa. S. aureus resistance to lincosamides can be acquired by modifications and/or mutations in the rRNA and rProteins. Here, we present the crystal structures of the large ribosomal subunit of S. aureus in complex with the lincosamides lincomycin and RB02, a novel semisynthetic derivative and discuss the biochemical aspects of the in vitro potency of various lincosamides. These results allow better understanding of the drugs selectivity as well as the importance of the various chemical moieties of the drug for binding and inhibition.

Ensemble- and Rigidity Theory-Based Perturbation Approach To Analyze Dynamic Allostery

Allostery describes the functional coupling between sites in biomolecules. Recently, the role of changes in protein dynamics for allosteric communication has been highlighted. A quantitative and predictive description of allostery is fundamental for understanding biological processes. Here, we integrate an ensemble-based perturbation approach with the analysis of biomolecular rigidity and flexibility to construct a model of dynamic allostery. Our model, by definition, excludes the possibility of conformational changes, evaluates static, not dynamic, properties of molecular systems, and describes allosteric effects due to ligand binding in terms of a novel free-energy measure. We validated our model on three distinct biomolecular systems: eglin c, protein tyrosine phosphatase 1B, and the lymphocyte function-associated antigen 1 domain. In all cases, it successfully identified key residues for signal transmission in very good agreement with the experiment. It correctly and quantitatively discriminated between positively or negatively cooperative effects for one of the systems. Our model should be a promising tool for the rational discovery of novel allosteric drugs.

Critical 23S rRNA interactions for macrolide-dependent ribosome stalling on the ErmCL nascent peptide chain

The nascent peptide exit tunnel has recently been identified as a functional region of ribosomes contributing to translation regulation and co-translational protein folding. Inducible expression of the erm resistance genes depends on ribosome stalling at specific codons of an upstream open reading frame in the presence of an exit tunnel-bound macrolide antibiotic. The molecular basis for this translation arrest is still not fully understood. Here, we used a nucleotide analog interference approach to unravel important functional groups on 23S rRNA residues in the ribosomal exit tunnel for ribosome stalling on the ErmC leader peptide. By replacing single nucleobase functional groups or even single atoms we were able to demonstrate the importance of A2062, A2503 and U2586 for drug-dependent ribosome stalling. Our data show that the universally conserved A2062 and A2503 are capable of forming a non-Watson–Crick base pair that is critical for sensing and transmitting the stalling signal from the exit tunnel back to the peptidyl transferase center of the ribosome. The nucleobases of A2062, A2503 as well as U2586 do not contribute significantly to the overall mechanism of protein biosynthesis, yet their elaborate role for co-translational monitoring of nascent peptide chains inside the exit tunnel can explain their evolutionary conservation.

Identification of potential allosteric communication pathways between functional sites of the bacterial ribosome by graph and elastic network models

BACKGROUND

Accumulated evidence indicates that bacterial ribosome employs allostery throughout its structure for protein synthesis. The nature of the allosteric communication between remote functional sites remains unclear, but the contact topology and dynamics of residues may play role in transmission of a perturbation to distant sites.

METHODS/RESULTS

We employ two computationally efficient approaches - graph and elastic network modeling to gain insights about the allosteric communication in ribosome. Using graph representation of the structure, we perform k-shortest pathways analysis between peptidyl transferase center-ribosomal tunnel, decoding center-peptidyl transferase center - previously reported functional sites having allosteric communication. Detailed analysis on intact structures points to common and alternative shortest pathways preferred by different states of translation. All shortest pathways capture drug target sites and allosterically important regions. Elastic network model further reveals that residues along all pathways have the ability of quickly establishing pair-wise communication and to help the propagation of a perturbation in long-ranges during functional motions of the complex.

CONCLUSIONS

Contact topology and inherent dynamics of ribosome configure potential communication pathways between functional sites in different translation states. Inter-subunit bridges B2a, B3 and P-tRNA come forward for their high potential in assisting allostery during translation. Especially B3 emerges as a potential druggable site.

GENERAL SIGNIFICANCE

This study indicates that the ribosome topology forms a basis for allosteric communication, which can be disrupted by novel drugs to kill drug-resistant bacteria. Our computationally efficient approach not only overlaps with experimental evidence on allosteric regulation in ribosome but also proposes new druggable sites.

Ligand-mediated and tertiary interactions cooperatively stabilize the P1 region in the guanine-sensing riboswitch

Riboswitches are genetic regulatory elements that control gene expression depending on ligand binding. The guanine-sensing riboswitch (Gsw) binds ligands at a three-way junction formed by paired regions P1, P2, and P3. Loops L2 and L3 cap the P2 and P3 helices and form tertiary interactions. Part of P1 belongs to the switching sequence dictating the fate of the mRNA. Previous studies revealed an intricate relationship between ligand binding and presence of the tertiary interactions, and between ligand binding and influence on the P1 region. However, no information is available on the interplay among these three main regions in Gsw. Here we show that stabilization of the L2-L3 region by tertiary interactions, and the ligand binding site by ligand binding, cooperatively influences the structural stability of terminal base pairs in the P1 region in the presence of Mg2+ ions. The results are based on molecular dynamics simulations with an aggregate simulation time of ~10 μs across multiple systems of the unbound state of the Gsw aptamer and a G37A/C61U mutant, and rigidity analyses. The results could explain why the three-way junction is a central structural element also in other riboswitches and how the cooperative effect could become contextual with respect to intracellular Mg2+ concentration. The results suggest that the transmission of allosteric information to P1 can be entropy-dominated.

Capturing RNA Folding Free Energy with Coarse-Grained Molecular Dynamics Simulations

We introduce a coarse-grained RNA model for molecular dynamics simulations, RACER (RnA CoarsE-gRained). RACER achieves accurate native structure prediction for a number of RNAs (average RMSD of 2.93 Å) and the sequence-specific variation of free energy is in excellent agreement with experimentally measured stabilities (R2 = 0.93). Using RACER, we identified hydrogen-bonding (or base pairing), base stacking, and electrostatic interactions as essential driving forces for RNA folding. Also, we found that separating pairing vs. stacking interactions allowed RACER to distinguish folded vs. unfolded states. In RACER, base pairing and stacking interactions each provide an approximate stability of 3-4 kcal/mol for an A-form helix. RACER was developed based on PDB structural statistics and experimental thermodynamic data. In contrast with previous work, RACER implements a novel effective vdW potential energy function, which led us to re-parameterize hydrogen bond and electrostatic potential energy functions. Further, RACER is validated and optimized using a simulated annealing protocol to generate potential energy vs. RMSD landscapes. Finally, RACER is tested using extensive equilibrium pulling simulations (0.86 ms total) on eleven RNA sequences (hairpins and duplexes).

Dual interaction of the Hsp70 J-protein cochaperone Zuotin with the 40S and 60S ribosomal subunits

Ribosome-associated J protein-Hsp70 chaperones promote nascent polypeptide folding and normal translational fidelity. Though known to span the ribosome subunits, understanding of J protein Zuo1 function is limited. New structural and crosslinking data allow more precise positioning of Saccharomyces cerevisiae Zuo1 near the 60S polypeptide exit site, pointing to interactions with ribosomal protein eL31 and 25S rRNA helix 24. The junction between the 60S-interacting and subunit-spanning helices is a hinge, positioning Zuo1 on the 40S, yet accommodating subunit rotation. Interaction between C-terminus of Zuo1 and 40S occurs via 18S rRNA expansion segment 12 (ES12) of helix 44, which originates at the decoding site. Deletions in either ES12 or C-terminus of Zuo1 alter stop codon readthrough and −1 frameshifting. Our study offers insight into how this cotranslational chaperone system may monitor decoding site activity and nascent polypeptide transit, thereby coordinating protein translation and folding.

In vivo quantification of the secretion rates of the hemolysin A Type I secretion system

Type 1 secretion systems (T1SS) of Gram-negative bacteria secrete a broad range of substrates into the extracellular space. Common to all substrates is a C-terminal secretion sequence and nonapeptide repeats in the C-terminal part that bind Ca²⁺ in the extracellular space, to trigger protein folding. Like all T1SS, the hemolysin A (HlyA) T1SS of Escherichia coli consists of an ABC transporter, a membrane fusion protein and an outer membrane protein allowing the one step translocation of the substrate across both membranes. Here, we analyzed the secretion rate of the HlyA T1SS. Our results demonstrate that the rate is independent of substrate-size and operates at a speed of approximately 16 amino acids per transporter per second. We also demonstrate that the rate is independent of the extracellular Ca²⁺ concentration raising the question of the driving force of substrate secretion by T1SS in general.

Predictive biophysical modeling and understanding of the dynamics of mRNA translation and its evolution

mRNA translation is the fundamental process of decoding the information encoded in mRNA molecules by the ribosome for the synthesis of proteins. The centrality of this process in various biomedical disciplines such as cell biology, evolution and biotechnology, encouraged the development of dozens of mathematical and computational models of translation in recent years. These models aimed at capturing various biophysical aspects of the process. The objective of this review is to survey these models, focusing on those based and/or validated on real large-scale genomic data. We consider aspects such as the complexity of the models, the biophysical aspects they regard and the predictions they may provide. Furthermore, we survey the central systems biology discoveries reported on their basis. This review demonstrates the fundamental advantages of employing computational biophysical translation models in general, and discusses the relative advantages of the different approaches and the challenges in the field.

Computational evidence that fast translation speed can increase the probability of cotranslational protein folding

Translation speed can affect the cotranslational folding of nascent peptide. Experimental observations have indicated that slowing down translation rates of codons can increase the probability of protein cotranslational folding. Recently, a kinetic modeling indicates that fast translation can also increase the probability of cotranslational protein folding by avoiding misfolded intermediates. We show that the villin headpiece subdomain HP35 is an ideal model to demonstrate this phenomenon. We studied cotranslational folding of HP35 with different fast translation speeds by all-atom molecular dynamics simulations and found that HP35 can fold along a well-defined pathway that passes the on-pathway intermediate but avoids the misfolded off-pathway intermediate in certain case. This greatly increases the probability of HP35 cotranslational folding and the approximate mean first passage time of folding into native state is about 1.67μs. Since we also considered the space-confined effect of the ribosomal exit tunnel on the cotranslational folding, our simulation results suggested alternative mechanism for the increasing of cotranslational folding probability by fast translation speed.

Multiple roles of the coding sequence 5' end in gene expression regulation

The codon composition of the coding sequence's (ORF) 5′ end first few dozen codons is known to be distinct to that of the rest of the ORF. Various explanations for the unusual codon distribution in this region have been proposed in recent years, and include, among others, novel regulatory mechanisms of translation initiation and elongation. However, due to the fact that many overlapping regulatory signals are suggested to be associated with this relatively short region, its research is challenging. Here, we review the currently known signals that appear in this region, the theories related to the way they regulate translation and affect the organismal fitness, and the debates they provoke.

Deletion of L4 domains reveals insights into the importance of ribosomal protein extensions in eukaryotic ribosome assembly

Numerous ribosomal proteins have a striking bipartite architecture: a globular body positioned on the ribosomal exterior and an internal loop buried deep into the rRNA core. In eukaryotes, a significant number of conserved r-proteins have evolved extra amino- or carboxy-terminal tail sequences, which thread across the solvent-exposed surface. The biological importance of these extended domains remains to be established. In this study, we have investigated the universally conserved internal loop and the eukaryote-specific extensions of yeast L4. We show that in contrast to findings with bacterial L4, deleting the internal loop of yeast L4 causes severely impaired growth and reduced levels of large ribosomal subunits. We further report that while depleting the entire L4 protein blocks early assembly steps in yeast, deletion of only its extended internal loop affects later steps in assembly, revealing a second role for L4 during ribosome biogenesis. Surprisingly, deletion of the entire eukaryote-specific carboxy-terminal tail of L4 has no effect on viability, production of 60S subunits, or translation. These unexpected observations provide impetus to further investigate the functions of ribosomal protein extensions, especially eukaryote-specific examples, in ribosome assembly and function.

Sequence-dependent elongation dynamics on macrolide-bound ribosomes

The traditional view of macrolide antibiotics as plugs inside the ribosomal nascent peptide exit tunnel (NPET) has lately been challenged in favor of a more complex, heterogeneous mechanism, where drug-peptide interactions determine the fate of a translating ribosome. To investigate these highly dynamic processes, we applied single-molecule tracking of elongating ribosomes during inhibition of elongation by erythromycin of several nascent chains, including ErmCL and H-NS, which were shown to be, respectively, sensitive and resistant to erythromycin. Peptide sequence-specific changes were observed in translation elongation dynamics in the presence of a macrolide-obstructed NPET. Elongation rates were not severely inhibited in general by the presence of the drug; instead, stalls or pauses were observed as abrupt events. The dynamic pathways of nascent-chain-dependent elongation pausing in the presence of macrolides determine the fate of the translating ribosome stalling or readthrough.

Deregulation of translation due to post-transcriptional modification of rRNA explains why erm genes are inducible

A key mechanism of bacterial resistance to macrolide antibiotics is the dimethylation of a nucleotide in the large ribosomal subunit by erythromycin resistance methyltransferases. The majority of erm genes are expressed only when the antibiotic is present and the erythromycin resistance methyltransferase activity is critical for the survival of bacteria. Although these genes were among the first discovered inducible resistance genes, the molecular basis for their inducibility has remained unknown. Here we show that erythromycin resistance methyltransferase expression reduces cell fitness. Modification of the nucleotide in the ribosomal tunnel skews the cellular proteome by deregulating the expression of a set of proteins. We further demonstrate that aberrant translation of specific proteins results from abnormal interactions of the nascent peptide with the erythromycin resistance methyltransferase-modified ribosomal tunnel. Our findings provide a plausible explanation why erm genes have evolved to be inducible and underscore the importance of nascent peptide recognition by the ribosome for generating a balanced cellular proteome.

Interaction of nascent chains with the ribosomal tunnel proteins Rpl4, Rpl17, and Rpl39 of Saccharomyces cerevisiae

As translation proceeds, nascent polypeptides pass through an exit tunnel that traverses the large ribosomal subunit. Three ribosomal proteins, termed Rpl4, Rpl17, and Rpl39 expose domains to the interior of the exit tunnel of eukaryotic ribosomes. Here we generated ribosome-bound nascent chains in a homologous yeast translation system to analyze contacts between the tunnel proteins and nascent chains. As model proteins we employed Dap2, which contains a hydrophobic signal anchor (SA) segment, and the chimera Dap2α, in which the SA was replaced with a hydrophilic segment, with the propensity to form an α-helix. Employing a newly developed FLAG exposure assay, we find that the nascent SA segment but not the hydrophilic segment adopted a stable, α-helical structure within the tunnel when the most C-terminal SA residue was separated by 14 residues from the peptidyl transferase center. Using UV cross-linking, antibodies specifically recognizing Rpl17 or Rpl39, and a His6-tagged version of Rpl4, we established that all three tunnel proteins of yeast contact the SA, whereas only Rpl4 and Rpl39 also contact the hydrophilic segment. Consistent with the localization of the tunnel exposed domains of Rpl17 and Rpl39, the SA was in contact with Rpl17 in the middle region and with Rpl39 in the exit region of the tunnel. In contrast, Rpl4 was in contact with nascent chain residues throughout the ribosomal tunnel.

Efficient and robust analysis of biomacromolecular flexibility using ensembles of network topologies based on fuzzy noncovalent constraints

We describe an approach (ENT(FNC)) for performing rigidity analyses of biomacromolecules on ensembles of network topologies (ENT) generated from a single input structure. The ENT is based on fuzzy noncovalent constraints, which considers thermal fluctuations of biomacromolecules without actually sampling conformations. Definitions for fuzzy noncovalent constraints were derived from persistency data from molecular dynamics (MD) simulations. A very good agreement between local flexibility and rigidity characteristics from ENT(FNC) and MD simulations-generated ensembles is found. Regarding global characteristics, convincing results were obtained when relative thermostabilities of citrate synthase and lipase A structures were computed. The ENT(FNC) approach significantly improves the robustness of rigidity analyses, is highly efficient, and does not require a protein-specific parameterization. Its low computational demand makes it especially valuable for the analysis of large data sets, e.g., for data-driven protein engineering.

Conformational readout of RNA by small ligands

RNA molecules have highly versatile structures that can fold into myriad conformations, providing many potential pockets for binding small molecules. The increasing number of available RNA structures, in complex with proteins, small ligands and in free form, enables the design of new therapeutically useful RNA-binding ligands. Here we studied RNA ligand complexes from 10 RNA groups extracted from the protein data bank (PDB), including adaptive and non-adaptive complexes. We analyzed the chemical, physical, structural and conformational properties of binding pockets around the ligand. Comparing the properties of ligand-binding pockets to the properties of computed pockets extracted from all available RNA structures and RNA-protein interfaces, revealed that ligand-binding pockets, mainly the adaptive pockets, are characterized by unique properties, specifically enriched in rare conformations of the nucleobase and the sugar pucker. Further, we demonstrate that nucleotides possessing the rare conformations are preferentially involved in direct interactions with the ligand. Overall, based on our comprehensive analysis of RNA-ligand complexes, we suggest that the unique conformations adopted by RNA nucleotides play an important role in RNA recognition by small ligands. We term the recognition of a binding site by a ligand via the unique RNA conformations "RNA conformational readout." We propose that "conformational readout" is a general way by which RNA binding pockets are recognized and selected from an ensemble of different RNA states.

Constraint Network Analysis (CNA): a Python software package for efficiently linking biomacromolecular structure, flexibility, (thermo-)stability, and function

For deriving maximal advantage from information on biomacromolecular flexibility and rigidity, results from rigidity analyses must be linked to biologically relevant characteristics of a structure. Here, we describe the Python-based software package Constraint Network Analysis (CNA) developed for this task. CNA functions as a front- and backend to the graph-based rigidity analysis software FIRST. CNA goes beyond the mere identification of flexible and rigid regions in a biomacromolecule in that it (I) provides a refined modeling of thermal unfolding simulations that also considers the temperature-dependence of hydrophobic tethers, (II) allows performing rigidity analyses on ensembles of network topologies, either generated from structural ensembles or by using the concept of fuzzy noncovalent constraints, and (III) computes a set of global and local indices for quantifying biomacromolecular stability. This leads to more robust results from rigidity analyses and extends the application domain of rigidity analyses in that phase transition points ("melting points") and unfolding nuclei ("structural weak spots") are determined automatically. Furthermore, CNA robustly handles small-molecule ligands in general. Such advancements are important for applying rigidity analysis to data-driven protein engineering and for estimating the influence of ligand molecules on biomacromolecular stability. CNA maintains the efficiency of FIRST such that the analysis of a single protein structure takes a few seconds for systems of several hundred residues on a single core. These features make CNA an interesting tool for linking biomacromolecular structure, flexibility, (thermo-)stability, and function. CNA is available from http://cpclab.uni-duesseldorf.de/software for nonprofit organizations.

Simulation study of the role of the ribosomal exit tunnel on protein folding

To investigate the role of the ribosomal exit tunnel on protein folding, we simulate the initial-stage folding behavior of the protein villin headpiece subdomain HP35 (PDB id: 1yrf) with and without prefolding in the exit tunnel by using an all-atom model and find that prefolding in the exit tunnel could effectively help the protein form native secondary structures. Furthermore, our results show that, after releasing from the exit tunnel, the prefolded chains may have a tendency to form more native contacts than those only in free space and this reduces the conformational space of sampling. Our results may provide an alternative way to explain the fast folding mechanism of proteins in vivo.

Structural basis of fluorescence quenching in caspase activatable-GFP

Apoptosis is critical for organismal homeostasis and a wide variety of diseases. Caspases are the ultimate executors of the apoptotic programmed cell death pathway. As caspases play such a central role in apoptosis, there is significant demand for technologies to monitor caspase function. We recently developed a caspase activatable-GFP (CA-GFP) reporter. CA-GFP is unique due to its "dark" state, where chromophore maturation of the GFP is inhibited by the presence of a C-terminal peptide. Here we show that chromophore maturation is prevented because CA-GFP does not fold into the robust β-barrel of GFP until the peptide has been cleaved by active caspase. Both CA-GFP and GFP₁₋₁₀ , a split form of GFP lacking the 11th strand, have similar secondary structure, different from mature GFP. A similar susceptibility to proteolytic digestion indicates that this shared structure is not the robust, fully formed GFP β-barrel. We have developed a model that suggests that as CA-GFP is translated in vivo it follows the same folding path as wild-type GFP; however, the presence of the appended peptide does not allow CA-GFP to form the barrel of the fully matured GFP. CA-GFP is therefore held in a "pro-folding" intermediate state until the peptide is released, allowing it to continue folding into the mature barrel geometry. This new understanding of the structural basis of the dark state of the CA-GFP reporter will enable manipulation of this mechanism in the development of reporter systems for any number of cellular processes involving proteases and potentially other enzymes.

Yeast polypeptide exit tunnel ribosomal proteins L17, L35 and L37 are necessary to recruit late-assembling factors required for 27SB pre-rRNA processing

Ribosome synthesis involves the coordinated folding and processing of pre-rRNAs with assembly of ribosomal proteins. In eukaryotes, these events are facilitated by trans-acting factors that propel ribosome maturation from the nucleolus to the cytoplasm. However, there is a gap in understanding how ribosomal proteins configure pre-ribosomes in vivo to enable processing to occur. Here, we have examined the role of adjacent yeast r-proteins L17, L35 and L37 in folding and processing of pre-rRNAs, and binding of other proteins within assembling ribosomes. These three essential ribosomal proteins, which surround the polypeptide exit tunnel, are required for 60S subunit formation as a consequence of their role in removal of the ITS2 spacer from 27SB pre-rRNA. L17-, L35- and L37-depleted cells exhibit turnover of aberrant pre-60S assembly intermediates. Although the structure of ITS2 does not appear to be grossly affected in their absence, these three ribosomal proteins are necessary for efficient recruitment of factors required for 27SB pre-rRNA processing, namely, Nsa2 and Nog2, which associate with pre-60S ribosomal particles containing 27SB pre-rRNAs. Altogether, these data support that L17, L35 and L37 are specifically required for a recruiting step immediately preceding removal of ITS2.

Global and local indices for characterizing biomolecular flexibility and rigidity

Understanding flexibility and rigidity characteristics of biomolecules is a prerequisite for understanding biomolecular structural stability and function. Computational methods have been implemented that directly characterize biomolecular flexibility and rigidity by constraint network analysis. For deriving maximal advantage from these analyses, their results need to be linked to biologically relevant characteristics of a structure. Such links are provided by global and local measures ("indices") of biomolecular flexibility and rigidity. To date, more than 14 indices are available with sometimes overlapping or only vague definitions. We present concise definitions of these indices, analyze the relation between, and the scope and limitations of them, and compare their informative value. For this, we probe the structural stability of the calcium binding protein α-lactalbumin as a showcase, both in the "ground state" and after perturbing the system by changing the network topology. In addition, we introduce three indices for the first time that extend the application domain of flexibility and rigidity analyses. The results allow us to provide guidelines for future studies suggesting which of these indices could best be used for analyzing, understanding, and quantifying structural features that are important for biomolecular stability and function. Finally, we make suggestions for proper index notations in future studies to prevent the misinterpretation and to facilitate the comparison of results obtained from flexibility and rigidity analyses.

Thermostabilizing mutations preferentially occur at structural weak spots with a high mutation ratio

We apply Constraint Network Analysis (CNA) to investigate the relationship between structural rigidity and thermostability of five citrate synthase (CS) structures over a temperature range from 37 °C to 100 °C. For the first time, we introduce an ensemble-based variant of CNA and model the temperature-dependence of hydrophobic interactions in the constraint network. A very good correlation between the predicted thermostabilities of CS and optimal growth temperatures of their source organisms (R²=0.88, p=0.017) is obtained, which validates that CNA is able to quantitatively discriminate between less and more thermostable proteins even within a series of orthologs. Structural weak spots on a less thermostable CS, predicted by CNA to be in the top 5% with respect to the frequency of occurrence over an ensemble, have a higher mutation ratio in a more thermostable CS than other sequence positions. Furthermore, highly ranked weak spots that are also highly conserved with respect to the amino acid type found at that sequence position are nevertheless found to be mutated in the more stable CS. As for mechanisms at an atomic level that lead to a reinforcement of weak spots in more stable CS, we observe that the thermophilic CS achieve a higher thermostability by better hydrogen bonding networks whereas hyperthermophilic CS incorporate more hydrophobic contacts to reach the same goal. Overall, these findings suggest that CNA can be applied as a pre-filter in data-driven protein engineering to focus on residues that are highly likely to improve thermostability upon mutation.

Flexibility analysis of biomacromolecules with application to computer-aided drug design

Flexibility characteristics of biomacromolecules can be efficiently determined down to the atomic level by a graph-theoretical technique as implemented in the FIRST (Floppy Inclusion and Rigid Substructure Topology) and ProFlex software packages. The method has been successfully applied to a series of protein and nucleic acid structures. Here, we describe practical guidelines for setting up and performing a flexibility analysis, discuss current bottlenecks of the approach, and provide sample applications as to how this technique can support computer-aided drug design approaches.

Kinetic analysis of ribosome-bound fluorescent proteins reveals an early, stable, cotranslational folding intermediate

Protein folding in cells reflects a delicate interplay between biophysical properties of the nascent polypeptide, the vectorial nature and rate of translation, molecular crowding, and cellular biosynthetic machinery. To better understand how this complex environment affects de novo folding pathways as they occur in the cell, we expressed β-barrel fluorescent proteins derived from GFP and RFP in an in vitro system that allows direct analysis of cotranslational folding intermediates. Quantitative analysis of ribosome-bound eCFP and mCherry fusion proteins revealed that productive folding exhibits a sharp threshold as the length of polypeptide from the C terminus to the ribosome peptidyltransferase center is increased. Fluorescence spectroscopy, urea denaturation, and limited protease digestion confirmed that sequestration of only 10-15 C-terminal residues within the ribosome exit tunnel effectively prevents stable barrel formation, whereas folding occurs unimpeded when the C terminus is extended beyond the ribosome exit site. Nascent FPs with 10 of the 11 β-strands outside the ribosome exit tunnel acquire a non-native conformation that is remarkably stable in diverse environments. Upon ribosome release, these structural intermediates fold efficiently with kinetics that are unaffected by the cytosolic crowding or cellular chaperones. Our results indicate that during synthesis, fluorescent protein folding is initiated cotranslationally via rapid formation of a highly stable, on-pathway structural intermediate and that the rate-limiting step of folding involves autonomous incorporation of the 11th β-strand into the mature barrel structure.

Role of antibiotic ligand in nascent peptide-dependent ribosome stalling

Specific nascent peptides in the ribosome exit tunnel can elicit translation arrest. Such ribosome stalling is used for regulation of expression of some bacterial and eukaryotic genes. The stalling is sensitive to additional cellular cues, most commonly the binding of specific small-molecular-weight cofactors to the ribosome. The role of cofactors in programmed translation arrest is unknown. By analyzing nascent peptide- and antibiotic-dependent ribosome stalling that controls inducible expression of antibiotic resistance genes in bacteria, we have found that the antibiotic is directly recognized as a part of the translation modulating signal. Even minute structural alterations preclude it from assisting in ribosome stalling, indicating the importance of precise molecular interactions of the drug with the ribosome. One of the sensors that monitor the structure of the antibiotic is the 23S rRNA residue C2610, whose mutation reduces the efficiency of nascent peptide- and antibiotic-dependent ribosome stalling. These findings establish a new paradigm of the role of the cofactor in programmed translation arrest in which a small molecule is recognized along with specific nascent peptide sequences as a composite structure that provokes arrest of translation. A similar mechanism could be used by the ribosome to sense a variety of cellular metabolites.

Computational approaches to RNA structure prediction, analysis, and design

RNA molecules are important cellular components involved in many fundamental biological processes. Understanding the mechanisms behind their functions requires RNA tertiary structure knowledge. Although modeling approaches for the study of RNA structures and dynamics lag behind efforts in protein folding, much progress has been achieved in the past two years. Here, we review recent advances in RNA folding algorithms, RNA tertiary motif discovery, applications of graph theory approaches to RNA structure and function, and in silico generation of RNA sequence pools for aptamer design. Advances within each area can be combined to impact many problems in RNA structure and function.

Nascent peptide in the ribosome exit tunnel affects functional properties of the A-site of the peptidyl transferase center

The ability to monitor the nascent peptide structure and to respond functionally to specific nascent peptide sequences is a fundamental property of the ribosome. An extreme manifestation of such response is nascent peptide-dependent ribosome stalling, involved in the regulation of gene expression. The molecular mechanisms of programmed translation arrest are unclear. By analyzing ribosome stalling at the regulatory cistron of the antibiotic resistance gene ermA, we uncovered a carefully orchestrated cooperation between the ribosomal exit tunnel and the A-site of the peptidyl transferase center (PTC) in halting translation. The presence of an inducing antibiotic and a specific nascent peptide in the exit tunnel abrogate the ability of the PTC to catalyze peptide bond formation with a particular subset of amino acids. The extent of the conferred A-site selectivity is modulated by the C-terminal segment of the nascent peptide, where the third-from-last residue plays a critical role.

SecM-stalled ribosomes adopt an altered geometry at the peptidyl transferase center

A structure of a ribosome stalled during translation of the SecM peptide provides insight into the mechanism by which the large subunit active site is inactivated.

As nascent polypeptide chains are synthesized, they pass through a tunnel in the large ribosomal subunit. Interaction between specific nascent chains and the ribosomal tunnel is used to induce translational stalling for the regulation of gene expression. One well-characterized example is the Escherichia coli SecM (secretion monitor) gene product, which induces stalling to up-regulate translation initiation of the downstream secA gene, which is needed for protein export. Although many of the key components of SecM and the ribosomal tunnel have been identified, understanding of the mechanism by which the peptidyl transferase center of the ribosome is inactivated has been lacking. Here we present a cryo-electron microscopy reconstruction of a SecM-stalled ribosome nascent chain complex at 5.6 Å. While no cascade of rRNA conformational changes is evident, this structure reveals the direct interaction between critical residues of SecM and the ribosomal tunnel. Moreover, a shift in the position of the tRNA–nascent peptide linkage of the SecM-tRNA provides a rationale for peptidyl transferase center silencing, conditional on the simultaneous presence of a Pro-tRNA^Pro in the ribosomal A-site. These results suggest a distinct allosteric mechanism of regulating translational elongation by the SecM stalling peptide.

In all cells, ribosomes perform the job of making proteins. As the proteins are synthesized they pass through a tunnel in the ribosome, and some growing proteins interact with the tunnel, leading to stalling of protein synthesis. Here, we used cryo-electron microscopy to determine the structure of a ribosome stalled during the translation of the Escherichia coli secretion monitor (SecM) polypeptide chain. The structure reveals the path of the SecM peptide through the tunnel as well as the sites of interaction with the tunnel components. Interestingly, the structure shows a shift in the position of the transfer RNA (tRNA) to which the growing SecM polypeptide chain is attached. Since peptide bond formation during protein synthesis requires precise placement of the substrates, namely, the peptidyl-tRNA and the incoming amino acyl-tRNA, it is proposed that this shift in the SecM-tRNA explains why peptide bond formation cannot occur and translation stalls.

Protein folding at the exit tunnel

Over five decades of research have yielded a large body of information on how purified proteins attain their native state when refolded in the test tube, starting from a chemically or thermally denatured state. Nevertheless, we still know little about how proteins fold and unfold in their natural biological habitat: the living cell. Indeed, a variety of cellular components, including molecular chaperones, the ribosome, and crowding of the intracellular medium, modulate folding mechanisms in physiologically relevant environments. This review focuses on the current state of knowledge in protein folding in the cell with emphasis on the early stage of a protein's life, as the nascent polypeptide traverses and emerges from the ribosomal tunnel. Given the vectorial nature of ribosome-assisted translation, the transient degree of chain elongation becomes a relevant variable expected to affect nascent protein foldability, aggregation propensity and extent of interaction with chaperones and the ribosome.

Transient tertiary structure formation within the ribosome exit port

The exit tunnel of the ribosome is commonly considered to be sufficiently narrow that co-translational folding can begin only when specific segments of nascent chains are fully extruded from the tunnel. Here we show, on the basis of molecular simulations and comparison with experiment, that the long-range contacts essential for initiating protein folding can form within a nascent chain when it reaches the last 20 Å of the exit tunnel. We further show that, in this "exit port", a significant proportion of native and non-native tertiary structure can form without steric overlap with the ribosome itself, and provide a library of structural elements that our simulations predict can form in the exit tunnel and is amenable to experimental testing. Our results show that these elements of folded tertiary structure form only transiently and are at their midpoints of stability at the boundary region between the inside and the outside of the tunnel. These findings provide a framework for interpreting a range of recent experimental studies of ribosome nascent chain complexes and for understanding key aspects of the nature of co-translational folding.

Recognition of the regulatory nascent chain TnaC by the ribosome

Regulatory nascent chains interact with the ribosomal exit tunnel and modulate their own translation. To characterize nascent chain recognition by the ribosome at the atomic level, extensive molecular dynamics simulations of TnaC, the leader peptide of the tryptophanase operon, inside the exit tunnel were performed for an aggregate time of 2.1 mus. The simulations, complemented by quantum chemistry calculations, suggest that the critical TnaC residue W12 is recognized by the ribosome via a cation-pi interaction, whereas TnaC's D16 forms salt bridges with ribosomal proteins. The simulations also show that TnaC-mediated translational arrest does not involve a swinging of ribosomal protein L22, as previously proposed. Furthermore, bioinformatic analyses and simulations suggest nascent chain elements that may prevent translational arrest in various organisms. Altogether, the current study unveils atomic-detail interactions that explain the role of elements of TnaC and the ribosome essential for translational arrest.

The key function of a conserved and modified rRNA residue in the ribosomal response to the nascent peptide

The ribosome is able to monitor the structure of the nascent peptide and can stall in response to specific peptide sequences. Such programmed stalling is used for the regulation of gene expression. The molecular mechanisms of the nascent-peptide recognition and ribosome stalling are unknown. We identified the conserved and posttranscriptionally modified 23S rRNA nucleotide m(2)A2503 located at the entrance of the ribosome exit tunnel as a key component of the ribosomal response mechanism. A2503 mutations abolish nascent-peptide-dependent stalling at the leader cistrons of several inducible antibiotic resistance genes and at the secM regulatory gene. Remarkably, lack of the C2 methylation of A2503 significantly function induction of expression of the ermC gene, indicating that the functional role of posttranscriptional modification is to fine-tune ribosome-nascent peptide interactions. Structural and biochemical evidence suggest that m(2)A2503 may act in concert with the previously identified nascent-peptide sensor, A2062, in the ribosome exit tunnel to relay the stalling signal to the peptidyl transferase centre.

Thermodynamic stability of polypeptides folding within modeled ribosomal exit tunnel: a density functional study

The mechanism of polypeptide folding, especially for the formation of tertiary structures, within the ribosomal exit tunnel, remains one of the most important unsolved problems in biophysical chemistry and molecular biology. In this work, we use a density functional theory (DFT) to explore the polypeptide folding within a modified nanopore, which mimics the confined environment of ribosomal exit tunnel. Results indicate that too long polypeptides (N>100 cannot fold into a helix state within the nanopore, and the helix polypeptides favor folding into a negative coiled coil rather than a positive one, because the negative coiled coil has a lower grand potential than the positive one, and the polypeptide folding into the negative coiled coil therefore needs less driving force than the positive one. To fold into the positive coiled coil, the helix polypeptides must have a small minor radius or a short chain length, which provides helpful insights into the design of nanodevices for manipulating the positive coiled coil. In the presence of attractive interaction, helices need more driving force to fold into coiled coil. Importantly, we have also proposed a scaling relation to understand the folding behavior. The scaling relation gives a good estimate for the computational results, and provides a reasonable explanation for the folding behavior. In summary, it is expected that the proposed DFT approach and the scaling relation provide alternative means for the investigation of polypeptide folding in confined environment, and these impressive results could give useful insights into nascent polypeptide folding.

Structural signatures of antibiotic binding sites on the ribosome

The ribosome represents a major target for antibacterial drugs. Being a complex molecular machine, it offers many potential sites for functional interference. The high-resolution structures of ribosome in complex with various antibiotics provide a unique data set for understanding the universal features of drug-binding pockets on the ribosome. In this work, we have analyzed the structural and evolutionary properties of 65 antibiotic binding sites (ABSs) in the ribosome. We compared these sites to similar-size computed pockets extracted from the small and large ribosomal subunits. Based on this analysis, we defined properties of the known drug-binding sites, which constitute the signature of a 'druggable' site. The most noticeable properties of the ABSs are prevalence of non-paired bases, a strong bias in favor of unusual syn conformation of the RNA bases and an unusual sugar pucker. We propose that despite the different geometric and chemical properties of diverse antibiotics, their binding sites tend to have common attributes that possibly reflect the potency of the pocket for binding small molecules. Finally, we utilized the ensemble of properties to derive a druggability index, which can be used in conjunction with site functionality information to identify new drug-binding sites on the ribosome.

Molecular recognition of RNA: challenges for modelling interactions and plasticity

There is growing interest in molecular recognition processes of RNA because of RNA's widespread involvement in biological processes. Computational approaches are increasingly used for analysing and predicting binding to RNA, fuelled by encouraging progress in developing simulation, free energy and docking methods for nucleic acids. These developments take into account challenges regarding the energetics of RNA-ligand binding, RNA plasticity, and the presence of water molecules and ions in the binding interface. Accordingly, we will detail advances in force field and scoring function development for molecular dynamics (MD) simulations, free energy computations and docking calculations of nucleic acid complexes. Furthermore, we present methods that can detect moving parts within RNA structures based on graph-theoretical approaches or normal mode analysis (NMA). As an example of the successful use of these developments, we will discuss recent structure-based drug design approaches that focus on the bacterial ribosomal A-site RNA as a drug target.

Divergent stalling sequences sense and control cellular physiology

Recent studies have identified several amino acid sequences that interact with the ribosomal interior components and arrest their own elongation. Whereas stalling of the inducible class depends on specific low-molecular weight compounds, that of the intrinsic class is released when the nascent chain is transported across or inserted into the membrane. The stalled ribosome alters messenger RNA secondary structure and thereby contributes to regulation of the cis-located target gene expression at different levels. The stalling sequences are divergent but likely to utilize non-uniform nature of the peptide bond formation reactions and are recruited relatively recently to different biological systems, possibly including those to be identified in forthcoming studies.